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GROWTH OF DISTURBANCES IN A FLAME-GENERATED SHEAR REGION ^ 

By Perry L. Blackshear, Jr. 


SUMMARY 

The growth of transverse velocity disturbances in the shear 
region caused by a flame in a duct was emmined theoretically 
and experimentally. In the theoretical stability analysis^ a 
flow fleld arising from a flame in a duct was analyzed. The 
flow was neutrally stable to symmetric disturbances and un- 
stable to antisymmetric ones. As a disturbance of given fre- 
quency 7noves downstream, from the flameholder ^ a most-amplified 
frequency is found that is equal to flow velocity divided by duct 
width. Disturbances at frequencies greater or less than this 
are less amplified. The amount that an incoming disturbance 
is amplified depends on the flame slope. For the most-amplified 
frequency y the ratio of terminal to initial transverse disturbance 
velocity is equal to where the slope is roughly turbu- 
lent flame speed divided by flow velocity. This ratio of ter7ninal 

to initial transverse disturbance velocity depends on inlet velocity 
andfla7ne speed only in the stated manner and is essentially in- 
dependent of density ratio. 

The critical wavelength sepamting disturbances that grvw 
frvrn those that do 7iot grvw varies as a constant times flame 
width. The constant is 2]i for flames that are 7iarrow cornparrd 
with duct width and becomes infinite as the flame reaches the 
wall. The flow is unstable to disturbances having wavelemgths 
greater than critical. Thus, short-wavelength (i. e., high- 
frequency) disturbances achieve their terminal amplifications 
near the fiarneh older, while long-wavelength disturbances achieve 
theirs fariher downstrvarn. The mvst-amplified disturbance 
achieves its terminal amplification when the flame fills 0.60 of 
the duct. The results are related to rneasumble quantities in 
the form of (1) flarne -front displacement, (2) phase velocity, and 
{3) to a first approxirnation, the e ffect of disturbance frequency 
on flame-propagation velocity. 

In the experimental pari of the prvgrarn disturbances of 
various frequencies were imposed on a flame stabilized in a 
duct, and the eflects werv measured by shadow photogmphy and 
2 )hotomultiplier-probe surveys. Kfforis to verify elements of 
the theory are encouruging but not conclusive. Reasonable 
agreement is achieved in the comparison of flame -frvnt displace- 
ment and phase velocity. In the measurements of the effect of 
frequency of disturbance on flame speed, there appeared to be 
a most-eflfective frequency, which was twice the most amplified 
frequency given by theory. There is no direct evidence of 
transition due to the amplified disturbances. The flow was 


turbulent in all cases owing to transition in the flarneholder 
boundary layer, either on the flarneholder or near the point of 
separatiori. Therv was some evidence of an increased small- 
scale turbulence in the flame when long-wavelength disturbances 
werv present. In other words, the large-arnplitude, long- 
wavelength disturbances cause a sort of tr ansition in the already 
turbulent flow and amplify the smaller - scale disturbances . 
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INTRODUCTION 

The stability of the flow field arising from an anchored 
flame in a duct is of interest for two reasons: (1) Experimental 
investigations have shown tliat the turbulent flame sjieeds 
attained in a flame anchored in a duct depend more on the 
turbulence originating in the shear regions generated by the 
flarneholder and flame than on the turbulence in the incident 
gas stream; and (2) some transverse modes of burner reso- 
nance encountered have been accoinpanied by an increase in 
the rate of flame propagation. 

DEPENDENCE OF TURBULENT FLAME SPEED ON 
DISTURBANCES ORIGINATING IN SHEAR REGIONS 

The problem of describing the flow field that arises when 
a flame is anchored in a duct was first posed and solved by 
Scurlock (ref. 1). One result of this work was the relation 
between the fraction of mixture liurned and the local flame 
width. Scurlock used this result to measure turbulent flame 
speeds in an experiment originally designed to give the effect 
of inlet turbulence on flame speeds. He found that turbu- 
lence associated with two regions of shear flow masked, to a 
great extent, the effect of imposed stream turbulence. These 


> Supersedes NAC A Technical Note 3830 by Perry L. Blackshear, Jr , 1956. 
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two ivgioiis are the near wake of tlie flaineliolder ami the 
sheai- region associated with the fully developed flame in 
a duct. The latter shear region is caused by the pressure 
gradient acting on gases of different densities. Scuiiock 
calculated the velocity profile in this region. 

In flow ])ast cylinders without flames, eddies are encoun- 
tered. Scurlock found that the flame seated on a cylindrical 
Jlameholder prevented eddies from shedding from the flame- 
holder over the range of Reynolds numliers investigated (up 
to aiiproximately 15,000). Haddock (ref. 2) and Zukoski 
and Marble (ref. 3), in investigating the mechanism of flames 
seating on cylindrical rods, found that the boundary layer 
separating from either side of the rod underwent transition 
from laminar to turbulent flow. The transition point ap- 
proached the point of separation at a Reynolds number of 
about 10k References 2 and 3 show that transition in the 
flaim'holder boundary layer is one mechanism for producing 
turbulence in the sheai* i*egion near the flameholder. 

The role played by the second shear i*egion cited by Scur- 
lock as a source of turbulence has not been directly investi- 
gated. After obtaining experimental results similar to those 
of Scurlock, several investigators [e. g., refs. 4 (Wold) and 5] 
conclrnh'd that tin*bulence is generated in this region. The 
arguments leading to this conclusion are briefly indicated 
in tlie following discussion. 

Tlieories of turbulent flame propagation have been iiostu- 
lated for two types of tin*bulence: (1) turbulence with a scale 
that is small compared with laminar llame-front thickness, 
and (2) turbulence witli a scale that is large compared with 
laminar flame-front thickness (refs. G to 8). Experiments 
conducted in the abseiu'e of shear regions yielded flame 
speeds in good agreement with values predicted by the large- 
scale-turbulence theoiy (refs. 9 to 12). Experiments con- 
ducted in the jiresence of shear regions (refs. 4, 5, and 8) gave 
flame speeds higlier than those predicted by theory. The 
intei’pretation of this difference is the subject of a current 
controversy (refs. 4, 5, 8, 13, and 14). It is generally agreed 
tiiat, when an existing shear region is of the type associated 
with an anchored flame in a duct, the higher flame speed is 
due in part to turbulence amplified by or originating in this 
shear region. 

Rayleigh (ref. 15) observed that some jets and diffusion 
flames are sensitive to certain imposed transverse disturb- 
ances. Fi’oman inviscid stability analysis, he found a 
velocity profile unstable to antisymmetric disturbances but 
neutrally stable to symmetric wavy disturbances. This 
result agreed with his observations on sensitive flames and 
jets. Some of the limitations of Rayleiglds woi*k are con- 
sidered in the theoretical section of tliis report. The analyt- 
ical juethod used by Rayleigh is ai)plied in the present ex- 
amination of the instability of flame-zone flow to imposed 
disturbances. 

EFFECT OF COMBUSTOR OSCILLATIONS IN INCREASING 
RATE OF FLAME PROPAGATION 

A great variety of combustion-chamber resonances have 
been observed in jet-engine combustors (refs. IG to 19). 
These i*esonances most often involve some mode of acoustic 
oscillation in the combustion chamber. Reference IG dis- 
cusses modes that involve longitudinal gas and flame motion 


and reports a deterioration in combustor performance when 
these oscillations occur. Keference 17 discusses a type of 
resonance identified as a transverse mode of oscillation in 
the combustion chamber. Accompanying this type of oscil- 
lation was an apparent increase in combustion efficiency 
(ref. 17) and in the rate of flame spreading (ref. 18). This 
type of oscillation, with an acoustic wavelength controlled 
by duct dimensions, could possibly excite a wavy disturbance 
that would subsequently be amjilified by the flame-generated 
shear i*egion. 

Unfortunately, these transverse modes (called screech) 
also cause the destruction of the burner walls when the 
screech becomes severe enough (ref. 18). Reference 19 
shows pictures of the flame near the flameholder when a 
combustor was resonating in a transverse mode. A strong 
interaction was apparent in the form of flame wrinkles form- 
ing at tlie flameholder and growing in amplitude as the^^ 
moved downstream. It seems apparent that at least one 
of the ways the transverse oscillation increases the rate of 
flame spreading is by wrinkling and hence extending the 
interface between tlie burned and unburned gases. 

If the oscillation that disturbs the flame at the flameholder 
is su])])lied locally (by some suitable exciting mechanism) 
and not by the sometimes destructive combustor resonance, 
some of the benefits of screech can be realized without its 
penalties. The gains in flame speed should be greatest for 
those frequencies of excitation that are amplified b}^ the 
flame. It is hoped that, by examining the sensitivity of a 
flame in a duct to imposed disturbances, an excitation 
frequenc^^ may be found that gives the optimum increase in 
flame speed. 

This report is divided into three parts — theory, experi- 
ment, and discussion. Separate introductions are given for 
the first two parts describing the purpose of each. The dis- 
cussion weighs the results of the other parts and includes 
some extrapolations. 

SYMBOLS 

A, B, C, 1) constants of integration 
b width of shear region in Rayleiglf s problem 

b' width of plateau in _/ \_ velocity profile 

c complex phase velocity, ft/sec 

d flameholder width 

E voltage 

F fraction burned 

f frequency, cps 

G gain; 

G' gain;U = <(f?')'/''^ 

h flame-front displacement from mean, ft 

/q duct half-width, ft 

Ju flame half-width, ft 

K\ constant, dF/[d(r/Ai)] 

Ko constant, d/q/dr, for an assumed flame-spreading 

mechanism 
])i*oport ionality factor 
k constant 

L Eulerian scale of turbulence, ft 

I length, ft 

If length of flame measured from wake, ft 

Iw length of flameholder wake, ft 
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F 

mean pressure, Ib/sq ft 

V 

pressure, Ib/sq ft 

V' 

])ci'turhation pressure, Ib/sq ft 

Re 

Reynolds number 

S 

flame surface area 

t 

time, sec 

u 

flow velocity in ^-direction, ft/sec 

l\ 

flow velocity in a^-direction in cold gas, ft /sec 

u* 

flow velocity in x-direction at apex of hot-gas 


velocity profile, ft /sec 

u' 

disturl)ance velocity in x-direction, ft /sec 

V 

flame speed, ft/sec 

v’ 

disturbance velocity in ?/-direction, ft/sec 

X 

coordinate in direction of duct axis, ft 

y 

coordinate perpendicular to duct axis, ft 

a 

wave number of disturbance, ft“^ 

OLo 

wave number of disturbance based on initial 


velocity and frequency, /27r/[/o 

OLs 

wave number for neutral stability based on 


inviscid analysis, ft"^ 

§ 

hot-gas velocity gradient, (f ^* — f i)//^ 2 ? ^ 

Iv 

e-c-v 

7i 


72 

e~ah., 

yi+y'i 


yl-y\ 

X 

wavelength, ft 

Xs 

2ir la„ ft 


1+75 
1 7| 

p 

density, lb-sec7ft‘‘ 

Pi 

density of cold gas, lb-sec Vft^ 

p> 

density of hot gas, Ib-sec^ft*^ 

p' 

perturbation density, lb-sec7ff^ 

Subscripts: 

^•distribution of disturbance amplitude, ft7sec 

cr 

critical 

d 

disturbed 

j 

frequency 

i 

imaginary part 

L 

laminar 

max 

maximum 

0 

initial; plane of Hameholder 

IHI 

measured from flame photographs 

pm 

measured by pliotom.ultipher 

r 

T 

real part 
turbulent 

t 

total 

u 

undisturbed 

THEORY: 

STABILITY OF FLOW FOR A FLAME IN A 


DUCT 


INTRODUCTION 

STEADY-STATE FLOW FIELD OF A FLAME ANCHORED IN A DUCT 

Tlic problem of (Icseribiiig the flow field tbat exists when 
a flame is anchored in a duct was first posed and solved by 
Scurlock in 1948 (ref. 1). Subsequently, Tsicn (ref. 20) 
and others (ref. 21) elaborated on the problem. The follow- 


ins: table summarizes vSciuioek^s ami Tsien^s assumptions 
and results as they pertain to the present work: 


Scurlock 

Common 

Tsien 


Assumptions 


(1) Fluid is incompress- 

(1) Fluid is steady-state, 

(1) Two cases are treated. 

ible. 

two-dimensional, inviscid. 

compressible and incom- 


non -heat-con ducting, 

pressible. In the inc«cn- 


non-turbulent. 

pressible, Scurlock’s 
assumption (2) holds: in 
the compressible, density 
is allowed to vary in the 
direction of flow. 

(2) Density jumps from 

(2) Velocity, density, and 

(2) Density ratio is con- 

PI to p 2 across the flame 

pressure are uniform at 

stant across flame front 

front and is constant 

combustor inlet. 

at any station. 

in each region. 

(3) Static pressure is uni- 

(3) The velocity profile is 

form at any cross section. 

composed of four straight- 


(4) The flameholder is of 
negligible size and spans 
the combustor. 

(5) The flame front is of 
negligible thickness. 

(6) The chamber is of con- 

line segments. 


stant width. 


- 

Results 


(1) Velocity profile as a 

(1) Fraction burned as a 

(1) Peak hot-gas \elocity 

function of flame 

function of flame width. 

as a function of flame 

width: slightly con- 

cave to incoming gas). 


width. 

(2) Momentum pres- 

(2) Cold-gas velocity as a 

(2) Limiting inlet velocity 

sure drop as a fund ion 

function of flame width. 

beyond which complete 

of flame width (i. e., 


combustion cannot occur 

of fraction burned). 


as a function of density 
ratio (for compressible 
case). 


A stability anabasis is made herein of a parallel flow having 
a velocit,y profile approximating the one found by Scurlock 
by means of straight-line segments and a density distribution 
as taken by Scurlock (i. e., constant Pi/po). 

STABILITY PROBLEM 

The study of the origins of turliulent motion has received 
considerable attention in the past, particularly with respect 
to turbulence occurring in boundary layers. One way of 
treating the problem is to assume the existence of infinitesi- 
mal wavy disturbances in the flowing stream and to seek tlie 
conditions under which these disturbances are amplified 
(refs. 22 and 23). In the application of this method to the 
study of tlie stability of the laminar boundary layer, a 
major objective is to determine a minimum Reynolds number 
below which all disturbances decay and to determine the 
effect of velocity profile on this minimum Reynolds number. 

The amplification of wavy disturbances supposedly leads 
to the transition to turbulence under certain circumstances. 
The desirability of maintaining laminar flow on airfoils makes 
an understanding of transition an important practical 
problem. In 1941, Schubauer and Skramstad conducted an 
experiment in which the disturbances introduced into a 
laminar boundary layer were carefully controlled (ref. 24). 
Tliey found a relation between the wave numbers of ampli- 
fied, neutral, and damped disturbances and the Reynolds 
number that confirmed the predictions of the infinitesimal 
wavy disturbance theory in almost every important respect. 
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The ])ractical value of learning about the stability of a 
flow field created by a flame in a duet is in furnishing in- 
formation that will enalile eombustor designers to encourage 
the growtli of disturbances in this region, because turbulence 
generally aids flame spreading. The information needed is 
(1) a range of Reynolds numbers over which the flow is 
unstable and (2) tlie kind of disturbance that possesses the 
greatest rate of anijilification within this range. In the 
following qualitative discussion most flames of practical 
interest apjiear to be in a Reynolds number range where 
amplified disturbances always exist and where the wave- 
lengths of disturbances that will be amplified may beobtained 
when viscosity is negle(*ted. 

If the disturbance varies as v' = where the 

wave numl)er a is real and the phase velocity c is generally 
complex, then the stability of boundary-layer flows (or of 
flows with symmetric profiles in a two-dimensional channel) 
can be illustrated by a grajih of values for a for neutral 
stability (i. e., where the imaginary part of c vanishes) 
against Reynolds number Re (ref. 22). 

A profile without a point of inflection would possess the 
following plot of a against Re for Ci={): 



Sketch (a) 

This curve is called the neutral envelope for the pi’ofile in 
question. Inside the envelope, disturbances are amplified 
(Ct>0); outside, disturbances are dampcal (Ci<0). Proceed- 
ing to the right fi*om the minimum Reynolds number attained 
1)3’' the curve are two branches of the curve As 

Re becomes indefinite! v gieat, both curves asvmptoticallv^ 
approach a=0. 

For the profile with a point of inflection, a similar relation 
a(c,Re) is obtained : 



Sketch (b) 


The significant difference between this plot for a profde with 
an inflection point and sketch (a) for a profile without an 
inflection point is that the asvmptotic values of the two 
branches approach different limits for large values of Re. 
The upper branch approaches a = otsj the wave number for 
neutral stability obtained bv neglecting viscositv, and the 
lower branch approaches a = 0. The Revnolds number range 
over which the flow with an inflection point is unstable 
extends to indefinitelv large Reynolds numbers. The flow 
generated b}^ a flame in a duct has, in effect, an inflection 
point. 

it is (lesiral)le either bv example or direct calculation to 
determine the order of magnitude of the Re 3 uiolds number 
range over which the neutral a departs from a.,; that is, 
where viscositv has a dominating influence. In two examples 
in the literature, Pretsch (ref. 25) and Lessen (ref. 26) 
examined a flow with a flex for stabilit 3 ^ In both instances 
the minimum Reynolds number for Cf=() was of the order of 
10; the maximum Re, where a for Ci=0 differed from as cal- 
culated from inviscid theorv bv more than 2 percent, was of 
the order of 10^ The profile of reference 25 is for a boundary 
laver with pressure rise; that of reference 26 is for tlie free 
boundarv h^ver between streams of differing velocities. 

If these values are taken to suggest the order of ReAUiolds 
numl)er where the stabilitv of a profile with a flex is sensitive 
to the effects of viscositv, some idea of flame-zone sizes can 
be estimated for the range of velocities encountered in jet 
engines (i. e., 50 to 500 ft/sec). For example, at 50 feet per 
second, the range of flame width where viscous forces are 
important is 1/8000 to 1/80 inch. At 500 feet per second, 
the range is 1/80,000 to 1/300 inch. Most of the flame zones 
of practical interest are much wider than these dimensions 
and mav be considered well within the range of Revnolds 
numbers where disturbances are amplified and where 
viscositv efiects are small. Therefore, the following inviscid 
stabilit\^ analysis ma 3 ^ give meaningful information on the 
value of tlie wave numbers a where amplification can occur. 

Lin (ref. 22) and Heisenberg (ref. 27) have criticized the 
use of broken-line profiles because there is a separate solution 
for c for each corner of the profile. These objections might 
be met bv the following observations: 

(1 ) The stead v-state flow field solved by Scurlock possesses 
corners (places where the velocitA^ gradient jumps) and can 
be very closelv approximated by straight-line segments 
without introducing new cornei*s into the ])rofile. The jumps 
in velocity gradient that occur at the corners result from the 
assumption that the flame front is of negligible thickness 
compared with the width of the flame zone. 

(2) Tlie following results show that the flame zone is 
unstable to disturbances having wavelengths 2.5 times the 
flame-zone width and greater. Thus, any rounding of the 
corners over a distance that is small compared with the flame- 
zone width should not affect the results; to the long- 
wavelength disturbances that are amplified, the cornei-s will 
still appear sharp. 
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In order to examine the effect of cornei*s in a velocity 
profile, three profiles of a similar shape are shown in the 
following sketch : 



Each of these profiles has been examined for instability, and 
a wavelength for neutral stability has been computed. The 
scale of the abscissa is chosen so that the value of the neutral 
wavelength is the same for the three profiles. Tlie two 
broken-line profiles are from Rayleigh (ref. 15); the contin- 
uous profile, from Savic (ref. 28). These solutions correspond 
to the antisymmetric disturbance. The flow in each case 
is inviscid. 

A similar comparison can be made between the velocity 
profile fl^=l + sinl(27r?///)— (7 t/ 2)] treated briefly by Lin 
(ref. 22) and a corresponding broken-line approximation. 
Lin’s result yields asl=w^j3 for the sine profile. Equation (37) 
in the following analysis was used to obtain the dimensions 
of a profile having the same asl. The following two profiles 
are given for flow between parallel walls, v' vanishing at the 
walls : 



Since the widths of the profiles composed of broken-line 
segments are so nearly equal to those with continuous veloc- 
ity gradients (having identical neutral wave numbers), it is 
believed that the broken-line profile used herein does not 
involve the introduction of extraneous solutions. 

INTERPRETING RESULTS OF STABILITY ANALYSIS 

The stability analysis gives an amplification rate aCt and 
a phase velocity Cr for disturbances as functions of the wave 
number a and the fraction of duct occupied by flame ho/hi for 
any one station along the axis of the duct. From tliis informa- 
tion a stability map is constructed, consisting of a plot of 
contours of constant aCi on a grid of ahi against hz/hi. This 
is a generalized map, with no assumption regarding the rate 
of change of fiaine width wdth distance. Lines of constant 
frequency are shown on the stability map calculated from the 
phase velocity Cr. From the map and an assumjition of the 


nature of flame-spreading rates, the total amplification of a 
disturbance of given frequency is calculated. 

Some aspects of this information are reduced to measurable 
quantities that can subsequently be compared with experi- 
mental results. These quantities include 

(1) Qualitative dependence of turbulent flame-spreading 
rate on disturbance frequency 

(2) Relation between critical wave number and flame 
width 

(3) Qualitative behavior of amplified disturbances as 
aflected by the steady-state condition (velocity, 
fuel-air ratio) 

One section of the analysis deals with the effect of a fiame- 
holder of finite size on the stead}^-state profile and hence on 
tlie stability analysis. 

STABILITY ANALYSIS 


STEADY-STATE SOLUTION 


The equations for incompressible flow derived by Tsien 
(ref. 20) relating t/i, f/^, t/*, pi/p 2 , ho/hiy and F are listed 
first; these relations are based on the assumption that the 
velocity profile can be approximated by straight-line seg- 
ments: 

Apex velocity: 


Uo 


V-s[C0-‘] 


( 1 ) 


Hot-gas velocity distribution: 


y. 

Ih2 

( 2 ) 


Flame width: 


^-1 

IT 

Fa 


Fraction burned: 




(3) 


(4) 


The hot-gas distribution (eq. (2)) and the relation between 
holhx and UijUo (eq. (3)) are used in interpreting the results 
of the stability anal^^sis. In the analysis the hot-gas velocity 
distribution (eq. (2)) is represented by 

U=U,+ Qi,~y)^ (y>0) 


where 


U=Ui-\-{h2-\-y)^ (?/<^ 0 ) 


Tsien’s analysis says nothing about the distribution of J^jhi 
with distance x. This must come from further assumptions 
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regarding flame slope, flame sjieed, or rale of change of F 
with X. For the present it is sufficient to say that the rate 
of change of hojhi witli x is small; that is, 

^ fh 

where X is the wavelengtli of a disturbance of interest. 

The stability of the profile assumed by Tsien, which closely 
approximates that calculated by Scurlock, is treated as a 
j)arallel-llow problem. 

FIKST-ORDEU DISTURBANCE PROBLEM 

Two excellent summary pajiers are available on the sta- 
bility of jiarallel flows (refs. 22 and 23). Xo attempt to 
review the field is made herein. The method employed here 
dates back to some of Rayleigh’s early work (ref. 15) with 
very slight modification. The development of the disturb- 
ance equation that follows is from reference 29. 

The flow is assumed parallel, inviscid, incompressible, 
plane, and non-heat-conducting. Gravity is neglected. 
The equations of motion are 


d.-“4 

^ dll , du , du\ 

v57+“S;+'al,) 

(5) 

/ 

by ^ ' 

( dr , dr , dr\ 

(6) 


The fluid is assumed incompressible (but the density is 
allowed to vary with ?/), so that 


dp dp dp ^ 

df dx ' dy 


(7) 


and the equation of continuity" is 


Let 


dx 


+^=() 


v=D'{x,y,f) ^ 

l>=P+p’{x,y,t) 

p=T{y)+p’{x,y,t) j 


(8) 


(5)) 


ia{U—c)p' + v' (12) 

— \~i(xu^ = Q ( 13 ) 

Eliniiiiating p', p', and u' from e((uatioiis (10) to (13) 
3 'i('kls an ocuiatioii for v' : 


(U-c) 


/dv , A ,d2r,idpr,„ 


c) 


dF 

dy 


F 



= 0 
(14) 


It is apparent that, though equations (7) and (12) arc 
useful in defining the tyjie of incompressible flow field as- 
sumed, they arc not required to obtain the disturbance 
equation (14) unless it becomes necessary to specify p . 

Defining 

F = — 


e(| nation (14) becomes 



d^' , 1 dp 
(l?/2 ”"p dy 



(hp_ (IT 
dy dy 



(15) 


In regions wliere p and dfV^ly constant, equation (15) 
reduces to 

"V=o (10) 

with the solution 

iP=A,e-^+A2e~<^^ (17) ‘ 


The disturbance velocities given by definition and equa- 
tion (13) are the real parts of 


ii^ = (^iiae^^ — Ayae 

F= — ia{Aie^^+A26~^^ 


ia(.x-ct) 


)gfo(x-cO 


} 


BOUNDARY CONDITIONS 


(18) 


The geometry of the flow field is shown in figure 1. Equa- 
tion (16) and the solutions (eqs. (18)) apply to each of the 
regions 1 to 4 considered separately. The establishment of 
the ])Oundarv-value problem consists in assigning suitable 


wliere primed quantities are of first order of smallness. 
Assume u' , F , and p to vary with x and t as 

where a is real and c, in general, is complex. 

By substituting equation (9) in equations (5) to (8) and 
retaining terms of the first order, equations (5) to (8) become 



-p[ia(U-c)v'\ = -^^'^^ (llj 


Region 1 | 

f 



h\ 


Flame 

P\ 

r 

1 




. w 


U-U^ 

2 

?2 1 



3 




Flame 






7^ 

4 

Pa'-P\ 


U^--Ux 


Figure 1. — Geometry of flow field. 
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bouiulaiy conditions at the walls and at the houndaries of 
regions 1 to 4. These boundary conditions are: 

(1) v' vanishes at y=±.h,x. 

(2) v' is continuous at interfaces. 

(3) j)' is continuous at interfaces (or, equivalently, the 
complete eq. (14) integrated across each interface is satisfied). 

Let the equation of tlie interface at y=hnt be 

y=hni+h{x,t) (19) 


(In the following derivation the subscripts m and n represent 
conditions evaluated in a region m or n.) Tlien, to the first 
order of approximation, at tlie interface. 


, bh.jjbh 


( 20 ) 


where U is the mean velocit}^ at the interface. In order for 
the pressure to be continuous at the interface, it is necessary 
that 

d/i d Xi 

where Xi is along the interface, or equivalently in the 
O’,?/ -system. 


dx ‘ dy dx dx dy dx 
where at the interface 

dy dh 

dx dx 


From the equation of motion, 


where 



u=U-\-u' 
v—v' >■ 


( 21 ) 


(5) 

( 6 ) 


(9) 


By vetaining terms of zero- and first-order approximations, 


^P_ ^p' 
br b.r 

bp_ bp' 
by~ by 



( 22 ) 


or bpjbx and bpiby are first-order terms. But, from equa- 
tion (20), it is clear that bh/bx is of the first order so that 
the products {bp/by){bh/bx^ are of the second order, and 
equation (21) liecomes 


bx bx 


or 







bu', .r 

b/ bx^^’"' by J 


(24) 


at the interface between the regions ?n and n. Eciuation 
(24) results from a necessary condition for the pressure to 
lie continuous at surfaces wliei-e either p or d^Vdy jumps. 

In order to show that this condition is ecjuivaleiit to re- 
quiring tliat the integrated equation (14) be satisfied at the 
interface, equation (14) is written 


d 

dy 



{U- 




br 



' — c)oPc' = () 


(25) 


Integrating equation (25) across the intei*face gives 




bv'n 

"" ^yJ 


which, with tlie insertion of the continuity equation, becomes 







Ikit u' varies as so that e([uation (26) can be recog- 

nized as another way of writing equation (24). 

The boundary conditions can be written 


v[=i) 

at 

y=l>i 

(27a) 

v[=V 2 

at 

y=hi 

(27b) 

djh_^2 
dx dx 

at 

y=h. 

(27c) 

t / 

V2 = Vz 

at 

y=() 

(27(1) 

II 

o'l.S^ 

at 

y=() 

(27e) 

v'i=v'i 

at 

y=-lh 

(27f) 

b]H_^'i 
bx bx 

at 

y=-h 

(27g) 

»;=() 

at 

y=—lh 

(27h) 


where the subscripts 1 to 4 denote ilow regions (fig. 1). 


SOLUTION OF THE BOUND ARY- VALUE PROBLEM 


The solutions of equation (16) for regions 1 to 4 (fig. 1) 
are written as follows: 


Region 1: 

Region 2: (p 2 =Ri^“^+R 2 ^““" ^ 

Region 3: 

Region 4: 


( 28 ) 
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Inserting equations (28) into equations (27) results in eiglit equations in tlie eight unknowns ^li, ^ 2 , • • • ^ 2 - The 
requirement that Ai, /I 2 , . . . Do have nontrivial solutions is 


y\ 

yl 


0 
— 1 


0 

—7? 


P2 P2 L L OiJ 


0 

0 

0 

0 

0 


0 

0 

0 

-1 


0 

0 

0 

—1 


(r*-e)+^ -[(t^^*-r)-^] (f^*-^)+^ 


yi 


0 

0 

0 

-1 

0 


-7i 


0 

0 

0 

—1 

0 

1 


r (r,-c)-^l7i -F (f/>-c)+^1 (L\-c)yl ^ (U,-c) 

L OiJ [_ aj P2 P2 


yi 


1 . 


=0 


(29) 


Where yi = e and yo = e By assigning = (this merely alters the real part of c) and operating on rows 1, 2, 

3, 0, 7, and 8, ecpiation (29) ean l)e reduced to the determinant 


c-i+;b, 

« P2 
1 


ho0-C + ^ 
a 




0 

-1 


0 

— 1 


-(M-c--) -(/(2/3-c--) L/3-C+- 

\ a/ \ a/ a 

0 ((Blc-c-i'W, f2lc+c-« 

\ P2 P2 Oi 


= 0 


(30) 


wliore- 


. 7|±7l 
^ il—y\ 


On oxpiuiding and simplifying, ocpiation (30) becomes 

(*■ S+0"1 (f S «+ S S «-S]=' 


(31) 


\vliere 


One mol eati be seen to l>e whollv real: 


t=l±Tl 

^ i-Ti 


,=?(t £!+£)' 

« \ P2 / 


(32) 


This can he identified as the root belonging to a symmetric disturbance by assigning z)' = z; 3=0 at ?/=0 in equation (27). 
In this case equation (30) reduces to 


Oi po po \ a/ 

1 1 
which yields equation (32) as before, ddie remaining two roots, 


= 0 


(33) 
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h2a-\- f — {^ah2— 1) ± '\/r /i'2«+ f 

P2 \ L P2 

■>] 



K«^2?-i) 

2| 


1 


(34) 


can similarly be identified as roots belon§>;ing to an antisymmetric disturbance. One way of doing this is to set at i/— /i3 
equal to v'^ at y= — h -2 or (from eq. (28)) 


/h+7^2Ti--eiT|-(72=0 

Using equation (35) in place of equations (27f), (27g), and (27h) gives for the determinantal equation 


(35) 


1 

yl 

0 

0 

0 

0 

1 

yl 

—1 


0 

0 

P2 

— cyl 

P2 

('-«>■ 

-(‘’+s)'*’ 

0 

0 

0 

0 

1 

1 

-1 


0 

0 

-c+^ 

a 

a 

c+- 

a 

a 

0 

0 

1 

t| 


— 1 


= 0 


(36) 


with roots of equation (34) for c as before. That the root for tlie symmetric disturbance c= (^/a)[f(pi/p 2 )+^]“* is wliolly real 
implies neutral stability. For the symmetric disturbance, then, the inviscid analysis does not yield information on wave 
numbers where amplification occurs, but merely suggests tliat the profile is relatively stable to symmetric disturbances and 
neutrally stable in the limit of infinite Keynolds number. 

The roots corresponding to the antisymmetric disturbances (eq. (34)) are complex when 


4^f + ^ (Wl'2“l)J 


(37) 


The complex root having the positive imaginary part corresponds to an amplified disturbance. The denominator of 
tlie right side of equation (34), 2[f(pj/p2)?+ 1], wlien written in full. 


\P2 


+ 1 


) 


VP2 1— 6“^“*' 

is positive for h- 2 <ihi-, tliat is, for the flame zone contained within the duct. Therefore, the amplified root of equation (34) is 




'-4(f^$+l)(a/l2^-D 

c=-S 

a 


> 


(38) 


Since this is an inviscid analysis, the damped root of equation (34) can have no meaning (ref. 22) and is conseiiuently ignored. 

RESULTS 

I'lie disturbance amplitude grows with time as where aCj is (from eep (38)) 


MI 


/ }l2a + f — (jioth'2 1 ) 
/ p2 

2 

^ all2^ 1 


f-s+1 

P2 


(39) 
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]() 


sinco 


0 = 


ho ho 


The quant it y uiuUm* (ho I’adical in equation (39) is iiide- 
])(‘ud(Mit of Jocal values. Tliis ([uantity is 

lilotted ill no-ur(' 2 as a function 2aho for the following 
(‘ondit ions: 

Pi i 7 • ho 

— =1 h ?n i-=0 

P2 h\—^oo 

-‘= 0 , 7 ('==0.25, 0..50, 0.75 

P2 Ih 

The first set of conditions corresponds to Rayleigh’s solution 
for the stability of an unconfined jet. The second set of 
conditions re])resents density ratios that are expected in a 
flame in a duct for three flame widths. The values of the 
parameter aCi{2holAU) in figure 2 go from 0 at 2aho=^0 
through a maximum, then to 0 again at a finite value of 
2aho. For a given value of 2Ao/df ^, this jiarameter behaves 
as the amplification rate aCj. 

It is clear from figure 2 that, for constant density, the 
flow is more unstable the smaller ho/hi becomes. For con- 
stant ho/hu the flow is more unstable the smaller pi/po be- 
comes. The influence of density ratio and wall proximity 
on the stability of the subject profile is in keeping with ex- 
])erimental observations of another type of ])rofile with a 
flex — namely, the wake of a cylindrical rod. The stabilizing 
effect of density ratio is illustrated in reference 1, which 
notes that a flame anchored in the wake of the rod inhibits 



Figure 2. — Amplification-rate parameter for _y^\_ profile. 


vortex shedding. The stabilizing effect of wall proximity is 
illustrated in reference 30, which shows that wall proximity 
increases the minimum Reynolds number for the occurrence 
of shed vortices. 

The stabilizing influence of the walls is further illustrated 
in figure 3. Figures 3(a) and (h) show the values of 2afi2 
for neutral stability and maximum instability plotted against 



Fraction of cross section filled by flame zone, ^ 3/^1 

(a) Wave-number parameter for neutral stability. 

(b) M^ave-mimber parameter for maximum instability. 

(c) Amj)lificat ion-rate parameter for maximum instability. 

Figure 3. — Wave-number parameter for neutral stability and maxi- 
mum instability, and amplification-rate parameter at maximum 
instability for profile. 
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A 2 // 0 . Figure 3(c) shows the amplification parameter 
aCi( 2 h 2 lAU) for maximum amplification rate plotted against 
ho/hi. Both the amplification rate and the range of the 
parameter 2 aJi 2 diminish as the llame fills the duct. 

CONSTRUCTION OF STABILITY MAP 

In considering the effect of disturbance amplification in 
the flame zone, local growth-rate contours as a function of 
wave numlier and ffame-width ratio l^/hi present nearly all 
the important information. Such a plot is called the stability 
map of a flame in a duct. The map is constructed from the 
data of figures 2 and 3 and from the steady-state relation of 
U*IUo and Ui/'LIo as functions of A 2 A 1 (1) and (3)) in 

the following manner. Given: {U\jUo){h 2 lhi) , {L'*IUo) (fhl 
hi), and [aCi{2h2lAU)](2ah2, A 2 M 1 ) for a given pi/p 2 - The 
quantity desired is [aCi(2hi/Uo)] (2a/q, Ih/hi). This quantity 
is obtained through suitable cross plots of figures 2 and 3 and 
equations (1) and (3) by setting 2ahi = 2ah2Kh2lhi) : 



The resulting stability map is shown in figure 4. Contours 
of constant dimensionless amplification rates aCi{2hilUo) are 
shown as functions of the dimensionless wave number 2ahi 
and dimensionless flame width A 2 /A 1 for pjpo of 7 and 5. 
These contours are in a region bounded l)y 2ahi= ^ , 0; 
hjh^ = i); and the line of 2ahi for neutral stability (i.e., the 
critical 2ahi as a function of h 2 /hi). 


PHASE VELOCITY 

A disturbance that distorts the flow and therefore distorts 
the flame front has a wavelength X=27r/a. Such a dis- 
turbance if introduced at the ffameholder travels along the 
flame front at the phase velocity c,. Tliis can be seen from 
the variation of the disturbance with x distance and time; 
that is, 

2;' ^ia^X-Crt) 

In onlci- to follow the amplification rate of a given dis- 
turbance, it is necessary to follow a disturbance of constant 
frequency. The relation between frequency and wavelength 
is 

(41) 

d'he quantity from equation (38) is given by 


Cr — Uj ^ 1 

U*-Ui ali2 


1 ) 


(42) 


for 

j^/^2a-t-f^ (^0:^2— 1) J ““4 ^ 

and by 


jr , l) + -\/ 

Cr — Ui 1 P2 \ 

{iochi— 


[-41 


1 {ah-A-\) 

U*-Ui ah-. 

I 

1 


1 


( 43 ) 
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for 


h20L-\- f 

P2 




( a ^ 2 ?— 1)>0 


A plot of this quantity against ah 2 is given for pjp.y = 7 at 
holhi = 0.25, 0.50, and 0.75 in figure 5(a). Also shown in 
figure 5(a) for eomparison is (c,— fq)/(f/^*— ?7i) for pJpo = 
5^ Ii 2 lhi = 0.25 and for pJpo=\, h 2 lhi = (). 

Xote that 


Ivn 

och2 — K) 


Cr — Ut 

U*-Ui 


=0 


From the steady-state solutions, 


lim (U*—Ui) = 0 

h.2/hi—*0 


Therefore, at tlie flameholder Cr=Ui=Uo- Then it is eon- 
venient to relate the local wave number to the wave number 
at the flameholder by 

ot=oto — (44) 

Cr 


where is the wave number at the flameholder, which is 
related to the frequency of the (listuri)ance/ and Uo by ao = 
12 ttIUo. By means of equations (42), (43), (44), (1), and 
(8), values of a{ao, /r>//?i) were calculated; these are shown 
in figure 4 as lines of constant frecpiency. These lines enable 
one to follow the local am])litication rate of a disturbance of 



(a) Pha.^o- velocity parameter for amplified and neutral disturbances 

as function of wave-number parameter. 

(b) Ratio of phase velocity to cold-gas velocity at three frequencies as 

function of h 2 lhi for density ratio of 7. 

Figure 5. — Phase-v^elocity i)arameters. 


given frequency from tlie flameliolder to some station down- 
stream and are so employed in a later section. 

An additional calculation for the phase velocity was made 
for later comparison with experimental data. In figure 5(b) 
is shown a plot of c,/Lq against for p,/p 2 of 7 and 2aohi of 
8, 12, and 16. These data enable direct comparison of 
measured values of phase velocity with values predicted from 
the stability analysis. 

DISTURBANCE GROWTH 

A qualitative picture of the nature of the disturbances 
])ropagating at phase velocity is shown in figure 6. The 
streamline distributions are given for a neutral disturbance 
in a reference plane moving at the phase velocity c,; in this 
frame of reference, the streamlines are no longer a function 
of time. Motion is relative to the ^Vyes'^ of the disturliance 
centers. In the hot flame core, the gas moves faster than c,; 
in the colder gas outside the flame, the gas moves slower 
than Cr. 

To picture what happens to an amplified disturbance, 
imagine the ?/-amplitude of the wrinkles in figure 6 to be 


Flow 



Figure 6. — Streamlines of neutral disturbance moving at phase 

velocity. 


growing from right to left. In addition, the flame width is 
increasing gradually, as is the phase velocity; and therefore 
the X- and ?/-spacing of the 'hwes'' is increasing in the direction 
of flow until a neutral situation is reached. At each j-station 
the disturbance has a local growth rate aCi that acts over a 
short time A/. Then (since x is related to time by dx=cgl/), 
at any time t after the disturbance originated at tlie flame- 
holder, the amplitude is 


v' 



otc^df 


(45) 


wlioro v'o is tlio aniplitiido of the distiirbaiicc at tlie flame- 
holder on the centerline of the jirofile.^ 

Up to this point, the independent variables have been the 
dimensionless wave number 2ah2 and flame width In 

order to evaluate equation (45), it is necessary to use the 
relation 




Cr 


(46) 


and some assumption concerning the variation of Ao/Ai with x. 
Two simple assum])tions are used that bracket the variation 
of h 2 lh\ with X found by experiment and facilitate the integra- 

2 The linking of the parallel-flow solutions by conserving along the centerline is the sim- 
plest but not necessarily the best of several available choices. 
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tion of equation (45). 
flame is linear: 


These assumptions are (1) that the 


I hi 


=K, 


h, 


(47) 


and (2) that the change 
linear : 


in fraction burned with distance is 


1 ^ 


=K, 


(48) 


A qualitative comparison of the way these assumptions 
describe a flame surface is given in tlie following sketcli : 



Inserting equations (46) and (47) in equation (45) results 
in 



h2lhx 

0 


2^1 Uo h\\!K^ 

Uo Cr h, \ 


(49) 


The term in parentheses is called the gain Gj which is inde- 
pendent of the value of the propagation rate K 2 . 

A similar derivation can be made for the alternative form 
of heat release; that is, for 


which 3 delds 



2 /»i Uo 
~UoTr 



(50) 


where the term in parentheses G' is independent of the value 

oiK,. 

The values of G and G' can be obtained by graphical 
integration. The}" give relative values of amplification of a 
disturbance of given frequency when integrated along a 
constant-frequency line in figure 4. Figure 7 gives the 
results of such an integration carried from A 2 A 1 of 0 to 1 for 
a range of 2oiohi and for the following conditions: G is given 
as a function of 2aohi for pi/p 2 of 5 and 7, and G' is given 
as a function of 2aolh for pi/p 2 of 7. It is apparent that, for 
the conditions examined, there is a most-amplified frequency 
464529—58 3 



corresponding to a value of 2aohi^^. Recall i\mij=2otohiUol 
2w2hi or the most-amplified frequency is approximately 
f^Uol^hi. This frequency does not depend upon a rate 
of flame spreading; it was obtained by assuming that a 
mode of flame spreading is specified. Approximately the 
same frequency obtains for the two modes of flame spreading 
specified. 

Not only is the most-amplified frequency in figure 7 at 
Pj/p 2 of 5 very near that for pjpo of 7, but also the numerical 
values of the gains are very nearly the same. This latter 
occurrence is due to two opposing influences: For a given 
f/*— ?7i, the profile having the lower value of pi/po has a 
higher amplification rate, and a disturbance of given fre- 
quency is amplified over a wider range of h 2 /lh. On the 
other hand, at any value of A 2 /A 1 , U*—Ui is smaller for 
the smaller pi/p 2 . 

The growth of the disturbances for fractional values of 
h 2 /hi (or F) was calculated by taking the integral in equation 
(49) or (50) to fractional values of A 2 A 1 (or F). The results 
are shown in figure 8, in which the gains G and G' are plotted 
against h 2 lhi or F for values of 2aohi from 2 to 18. For some 
values of /i 2 Ai"^ 5 , there are frequencies higher than those 
corresponding to 2ao^i«6 that locally exhibit a maximum 
gain. The fact that for values of /i2Ai<^0.5 there is a local 
amplification that depends on frequency is used later in 
comparing analysis with experiment. As h 2 lhi decreases, 
the frequency that is most amplified increases. 
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Fraction of cross section filled by flame zone, 

I 1 I I I I I [ I - I 1111 

0 .05 .1 .15 .2 .25 .3 .4 .5 .6 .7 .8.9 I 


Fraction burned, F 

Figure 8. — Gains as fuiictioii of hilh\ for family of initial wavo-nninlxn* 
parameters. Density ratio, 7. 


DISTURBANCE VELOCITY DISTRIBUTION IN i/-DIRECTION 

The way tlie disturbance velocity v' varies along the axis 
of the duct is sliown in figure 8. In reducing tlie results 
of this analysis to measurable quantities, a qualitative 
picture of the way the disturbance velocity v' varies across 
the duct is needed. Such a picture can most easily be 
obtained by considering a neutrally stable disturbance. 

From equations (18) and (28), 


4. 


T// jy±\ 1 


()<7/</r. 


(51) 


where V 2 ,y is the disturbance velocity at y in region 2, is 
the disturbance velocitv at the centerline, and yy = e~^^. In 
the rc'gion Ih^y obt aim'd from ecpia- 

tions (18), (27), and (28) in tlie form 


7i+7y 72 
-~7 i+72 7i/ 


(52) 


The ratio BijBo oan be obtaiiu'd from the equation rejire- 
sented b}^ the first row in the determinantal equation (30) 


or 


\ a P 2 / \ a p 2 / 


\ P2 y 


' ac 

(f-+7 


LV P2 ^ 



Recall that ^/ac in equation (30) is the same as 

Ui)/{c—Ui)], which for real values of c is given in 
equations (42) and (43). 

A series of calculations was made to obtain the representa- 
tive behavior of across the duct, and the results are 

shown in figure 9 for pjp 2 of 7. The critical ^ 2 //^ is taken as 
the value where for a given frequency the imaginary part of 
c becomes 0. Figure 9(a) shows the value of at the 

critical plotted against initial wave number 2a„/ii. For 
2aohi >8 (i. e., for (/i2/Ai)cr<0.4()) the value of v'^Jv'^^0.2. 
The I’atio decreases at 2a^Ai<8 (i. e., for (A2/Ai)cr>0.46) as 
the flame approaches the wall. 

In figure 9(b) the value is given for a fixed frequency 
corresponding to 2aohi = \2 following a disturbance along the 
duct from the flameholder (h 2 /hi = i)) to the station where 
the flame touches the wall. The dashed portion is intei- 
polatcd between ho/hi = i) (where and / 12 A 1 = (^ 2 ^ 1 ) cr 

(where the calculation starts). It is obvious from this curve 
that, if is measured instead of the apparent growth 
rate is not as rapid for /^ 2 //q<^(A 2 //q)cr for v'^, simply be- 
cause of the velocity distribution. Also, for A 2 //^iXA 2 //q)cr 
tliere is an ajiparent decay in owing to this distribution, 
despite the fact that the value of is neutrally stable. This 
latter point is employed in comparing the analyth'al with the 
experimental results. 

44ie significance of figure 9 is summarized as follows: 

(1) Although the inviscid treatment eliminates considera- 
tion of damping, it is clear that at the flame front the velocity 
amplitude diminishes at /? 2 /^i>(A 2 /Ai)cr because of the veloc- 
ity distribution. 

(2) In some instances, even though the flow is unstable 
and the disturbance velocity amplitude is increasing along 




(a) At critical ho/hi. 

(b) For wave-number parameter of 12. 

Figure 9. — Disturbance-velocity distributions. Density ratio, 7. 
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the center of the flame, the disturbance-velocity amplitude 
at the flame front may not necessarily be increasing. 

(3) The disturbance that is most amplified at the center 
of the profile (i.e., 2aA = h) is not necessarily the most 
amplified at the flame front, owing to the velocity distribu- 
tion; the frequency of the latter is somewhat greater than 
that of the former. 

relation of theoretical results to measurable quantities 

In the experimental phase, current instrumentation tecli- 
niques impose limitations on what can be measured. For 
example, the velocit}^ cannot be measured directl}^ by hot- 
wire techniques because the flame temperature is too high 
(«4000° K). The purpose of tliis section is to collect cer- 
tain results of the foregoing anal 3 ^sis that lend themselves 
to direct measurement. 

Flame slope. — In equations (49) and (50) the growth of a 
disturbance v' is found to be proportional to a calculable 
({uantit}^, gain, raised to a power \jK. The quantity K is 
alternatively /i 2 =d(A 2 //h)AK‘^Ai) or Ki = dFld{xlhi)^ depend- 
ing on the mode of flame propagation assumed. For both 
modes, K is smaller the more gradual the rate of flame 
spreading. To a first approximation, then, the amplification 
of a disturbance increases as the ratio of propagation velocity 
to flow velocity decreases. This trend can be verified experi- 
mentall3^ 

Phase velocity. — One quantity that can be measured 
immediately is the phase velocity Cr. A plot of CrjUi for tln-ee 
frequencies is given in figure 5(b). Thus, for a given flow 
setting, the absolute magnitudes of CrjUi can be compared 
and the variation with frequenc}^ and with /i 2 //q can be ex- 
amined. 

Flame-front displacement. — The transverse velocity v' 
acting at the flame surface causes a wavy’* wrinkle to 

appear in the surface. If a quantity h is defined as the 
distance from crest to valle.y of these wrinkles, this quantity 
equals twice the y amplitude of the streamlines of the flow 
at tlie flame front. If a neutral case is considered and the 
distribution of r' in the flame zone is approximated by 

— (l~^) 

tlie equations of the streamlines are given by 

— € cos ar± Vl“2^*+2fc cos ax 
Vc 


where 


Vc 


Vc 

K 

Vc 


{U*—c)ayc 

amplitude of disturbance velocity at ijc 
amplitude of disturbance velocity” at centerline 
coordinate where U=c 


Vo — Vc 


(U^—c)ayc 

T* dimensionless stream function. 


{U*—c)y, 


£udy 


At h2 ^ h2 , cr , S’F* varies between 0.96 and 0.99. Because 2 k 


can be taken arbitrarily” small, for 2/:’/(l — 24^*)<d, 


V 


— 1 — e COS ax 




Because y^yc, tlie plus sign is clioseii so that 


■^=1+^1— 2^*+ 


(v(l-2^*) *) 


cos ax 


<l-t-0.2+ 






Uo Uo 


U*-U, (0.2) 


U„ 


U„ 


cos aX 


' 1 . 2 - 


Uo 


U* 


■—c\v„ / 


COS aX 


But Vclv'o is a weak function of a and ^ 2 //^ and is approxi- 
mately” 0.2. 

Employ”ing the identity” 

f/* r /U* U^\f c-U, \ 

Uo Uo\Uo Uo)\ u^-uj 

and a value of {c—Ui)l{U*—Ui) from figure 5(a), 

U* c /U* Ud 


Uo Uo~\Uo Uo)^^^^ 


where the quantity {U*IUo) — {UilU„) is a function of p,/p 2 
and A 2 A 1 oiity. Thus, 




0 — 
^■^Uo 


'o Q.. 
0.85 aye 


■ cos ax 


at the flame front. The coefficient of cos ax is \ h , the flame- 
front wrinkle amplitude. It can be seen that h will heliavc 
((ualitatively like v' . Without going into further detail, then, 
some predictions can he made from the theoiw that may he 
verified by e.xperiment. First, h may be e.xpected to behave 
(lualitatively as shown in tlie following sketch: 



Sketch (f) 

The shape of the curve in sketch (f) can be used to locate 
the appro.ximate x-position of {h 2 /hi)cr- There is a jioint ol 
inflection where reaches a ma.ximum and a knee as 

grow. The knee position falls between the 
LJhi for maximum amplification and {hilhi)^- Thus, the 
data that can be sought are the absolute values of h 2 lfh 
where the knee occurs for a given frequency and the relative 
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positions of tlio knee for (lifforont frequencies. 

Flame propagation. —As stated in the introduetion, one of 
the reasons for pursuing the stability analysis is to find a 
disturbanee frecpieney that permits the greatest increase in 
flame spreading with the smallest expenditure of energy. In 
this section a model of flame propagation is constructed that 
permits extrapolating the behavior of v' from the staliility 
analysis to a measurahle behavior of flame propagation. 
The flame front is assumed to be a zone in which the density 
goes from pi to po in a distance that is short compared with 
the disturbance wavelength. This front is propagating into 
the gas of density pi at a rate Vt.u] is, the undisturbed 
turbulent flame speed. Owing to the normal ])ropagation 
of the front, an initially wavy surface becomes a cusped 
surface. Assuming sucli a surface composed of cylindrical 
arcs 


Unburned 



Burned 


Sketch (g) 


permits a calculation of the flame-front stretcliing due to a 
displacement h. If the undisturbed area X and tlie 
disturbed area Sa are approximately equal to arc lengtli, 





An equilibrium h is reached when 


m 




Vt, 

tilde. Equation (54) contains a nonlinear term 4/^VX^ 
The development so far as l)een a linear one. The linear 
development would be expected to furnish only a first 
approximation for amplitudes large enough so that 4h^/\^ 
is no longer negligible compared to unity. 

The sequence of events that leads to a prediction of Sa/S^ 
is as follows: 

(1) The flame is initially wrinkled by the velocity v' 
acting at the flame surface. 

(2) As the wrinkles grow, the flame shape shifts from a 
sine wave to a cusped surface, owing to normal propagation. 

(3) After a certain time an equilibrium is reached between 
the disturbance velocit}^ v' and the undisturbed propagation 
velocity Vt,u- 

The ratio of the disturl)ed surface to the undisturbed 
surface when equilibrium is reached is, for small disturb- 
ances, 


- j for small ampli- 
uJ 


Then Sa\Su behaves roughly as 








(56) 


Before this equilibrium value of Sa\Su is reached, the flame 
surface is given by (again for the neutral case) 


Vc 

But 2hj\ of equation (54) is given b}’ 


cos ad' 


2h 

X 



7T 


Vg 0.4 

AT 7T 


cos ax 



Figure 10. — Frequencies that sive maximum increase in local flame 
speed. Density ratio, 7. 
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Thus, in the transient case as in the equilibrium case (eq. 
(55)), SJSu beliaves as v\ Tlie quantity Sa/Su when meas- 
ured locally should behave like G in figure 8, if the effect of 
distribution is neglected. When the latter effect is taken 
into account, v' for fixed 2oLohi measured at h 2 achieves a 
maximum value somewhere between A 2 A 1 for maximum 
amplification and neutral stability. 

Thus, at any station h 2 lhi, there should be a frequency of 
disturbance that would cause the greatest increase in flame 
speed, and at this station the increase in flame speed should 
be greater for greater values of v' IV t,u- 

Figure 10 gives a ])lot of 2aohi against /i. 2 //q for maximum 
ani])lification and for neutral stability. For any value of 
/i, 2 //q, the frequency giving the greatest v'IVt,u (jd tbe flame 
surface) has a value of 2«y/q that falls between these curves. 
Tlie lag in time required for the flame to become wrinkled 
initially will move the curve for maximum Sa/Su to slightly 
higher 2aohi values. The curves given should be considered 
the lower limits of 2aohi for maximum Sa/Su, and they serve 
as a basis for comparison with the experimental data. 

KFFECT OF FLAMEHOLDER 

When a dame is anchored in a duct, some form of flame- 
holder must be used. The most common form of flame- 
holder studied in recent years is the cylinder with axis 
normal to the stream. The cross section of the cylinder has 
assumed many different shapes. Properties of the flow in 
the wake are common to all the shapes studied. In this 
section some of these properties are used to deduce the 
effect of a flameholder on the How field and on its stability 
far from the flameholder. 

Steady-state velocity distribution. -The properties of the 
flow held in the immediate wake of the flameholder as 
described in references 1 and 31 are shown in sketch (h): 

-Attached recirculation zone 



Two eddies are attached in the wake of the flameholder. 
The gas in this zone is hot; it acquires energy from the down- 
stream hot gas and delivers energy to the outside cold gas 
by diffusion. At some position on the surface of the recircu- 
lation zone, sufficient energy is acquired by the cold gas for 
ignition to occur and a flame to propagate into the cold-gas 
stream as shown in sketch (h). The pressure measured at 
the walls has an axial distribution that is shown qualitatively 
in sketch (h) (refs. 1 and 31). 

If the static pressure is assumed constant across tlie duct 
at any axial station, the effect of the first positive, then 
vanishing, then negative pressure gradient upon gases of 
unlike densities permits an estimation of the velocity pro- 
file downstream of the flameholder. The velocity profiles 


at stations X, Y, and Z of sketch (h) are shown in sketch (i): 



At station X the pressure gradient in the axial direction has 
been positive for a small distance and is still positive at X. 
This fixes the sense of the velocity difference between cold 
and hot gases on either side of the flame and the sense of 
the velocity gradient of the hot gas at the flame front. At 
station Y the gas has a jiositive-pressure-gradient history, 
again fixing the sense of the velocity difference, and at Y 
the pressure gradient vanishes, causing the velocity gradient 
at the flame front to vanish. At station Z the gas has a 
negative-pressure-gradient history from Y sufficient to cause 
the sense of velocity difference and gradient at the flame 
front as shown. The gas particle at the position of the 
velocity maximum passed through the flame at station Y. 

If the profile is approximated by straight-line segments 
wit h a corner at this velocity maximum as shown in sketch (j). 



some quantitative relation based on the continuity equation 
and experimental observations can be obtained. From con- 


tinuity. 

1 b'U*=VJ^^ 

2 P2 


and 

(57) 

where 

,U*+l\ p. 

0 ^ V pip 

^ P2 

(58) 


Ip axial length from station Y to station of interest down- 
stream of Y (sketch (h)) 
axial length of wake to station Y 
Vf mean flame speed over Ip 
mean flame speed over 

Combining equations (57) and (58) gives 

6 ' U^+U, VJ^ 

26” 2t7* Vplp ^ ^ 
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In order to obtain an order of magnitude of V I2h of interest 
in the stability analysis, the terms on the right side of equa- 
tion (59) are given the following values: (J7*+ C/i)/?7* = 1.5 
(from ref. 1 for ?'=0.25, Pi/p 2 = 6) ; T 1/(1 VI l) ~ 1/2.5; 

and (from ref. 1). Defining I as 


L 

t p 


~1 


equation (59) becomes 
b' 

26" 


1.2 




(60) 


At a station ljd^% (from eq. (60)), 



Farther downstream, at //c?=24, for example, 



The foregoing stability analysis is for a profile in which 6' = 0. 

From this section it can be seen that, near the flameholder 
(//^/=8), the hot-gas velocity distribution differs in shape 
from tlie _/\_ profile neglecting the flameholder. At this sta- 
tion, 60 percent of the flame-zone width has a flat distribu- 
tion. Kemote from the flameholder (//r/=24), 6 percent of 
the hot-gas width has a flat velocity distribution. The mag- 
nitude and direction of the error incurred in neglecting this 
flattening are estimated in the following section. 

Comparison of stability of_/\_and_/ \_ profiles. — Ampli- 
fication rates as a function of 2aho are compared in figure 2 
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Figure J 1. — Effect of velocity profile on amplification-rate parameter 
and on wave-number parameters for neutral stability and maximum 
instability. 


for several values of pi/p 2 and hilhx. The effect of increasing 
Pi/p 2 and h2lhi is stabilizing; that is, the magnitude of the 
amplification rates and the range of 2ah2 over which positive 
values exist are reduced. The plots of amplification rate 
against 2ahi are affine. It appears reasonable that, in exam- 
ining the effect of truncating the profile, an examination at 
a single value of pi/p 2 and would yield representative 

results. 

The conditions considered here are the ones formulated by 
Rayleigh in equations (24) and (25) (p. 387, ref. 15), which 
correspond to p\jp2= 1 and A 2 /^i = 6 (or hi = <x>). Calculations 
were made in which 6'/26 was taken as a parameter in such 
a way that the results 3 ueld values of aCi(26/At/) and 2ab at 
maximum instability and a critical value of 2ab where aCt 
vanishes. The results of these calculations are shown in 
figure 11. When 6' = 0, the results correspond to those of 
figure 2; the flow is unstalile to antisymmetric disturbances 
and neutralty stable to syunmetric disturbances. As b' j2b 
increases, the amplification rate and the critical values of 2ab 
decrease for the antisymimetric disturbances. The flow 
gradually becomes unstable to svmmetric disturbances as 
b' j2b increases. At 6726<C6.2, the flow is as unstable to 
symmetric as to antisvmmetric disturbances and possesses 
the same critical values of 2ccb. The values of b' j2b of interest 
here are of the order 0.6 to 0.06. The flame half-width is 
6'/2 + 6. In the following table the values of the important 
parameters aCi(2h2lAU), (2ah2)max cp find {2ah2)cr are com- 
pared for b'/2b = 0, 0.06, and 0.6 for antisymmetric dis- 
turbances: 


1/ 

2b 


(2cth2)maxei 

(20^2) er 

0 

0. 492 

2. 44 

3.69 

.on 

.508 

2. 49 

3. 71 

.() 

.688 

3. 28 

4. 57 


For a given 2li2lAUj the amplification rate is increased as 
the hot-gas velocity^ is jirogressivelv truncated. The wave- 
number range over which the flow is unstable is widened. 
Both of these eflects bespeak a greater instability' of flow. 
On the other hand, it can be easily' demonstrated that, for 
a given AU^ is greater w'hen the profile is truncated, so 
that the destabilizing effect of the flameholder is partially' 
but not completely attenuated. 

The difference between the parameters calculated at 
6'/26 = 0 and 0.6 shows that near the flameholder the amplifi- 
cation rates may' be of the order of 20 percent higher for 
some wave numbers than were those previously found. 
Because the flame zone must start with a finite width depend- 
ent upon the flameholder width, those frequencies having 
wave numbers 2 ad that exceed 2 ad sup- 

pressed. Here, JT is some proportionality factor less than 
unity. Certainly', neglect of the effect of a finite-size 
flameholder facilitates generalizing the results of the stability 
anah'sis. The investigations of this section indicate that, 
for flameholders small compared with duct width, such 
neglect appears warranted in a large part of the flame zone. 
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SUMMARY OF THEORETICAL RESULTS 

1 . The flow field of a flame in a duct is unstable to antisym- 
metric disturbances having wavelengths of the order of 2.5 
times the local flame width and greater. As A 2 A 1 approaches 
unity, this critical wavelength approaches 00 . 

2. The amplification rate depends directly upon the mag- 
nitude of the velocity difference (apex velocity minus cold-gas 
velocity) and inversely upon the flame width. Foi' fixed 
values of these quantities, the maximiun amplification rate 
decreases with increasing values of pi/p 2 and A 2 A 1 , going to 
zero for ^ 2 //^= 1 • 

When the history of a disturbance is followed through 
the course of an entire flame, a wave number for maximum 
disturbance growth is found. This maximum growth is for 
such a wave number at the flameholder that 2a,Ai = (), which 
corresponds to a frequency /^ f The actual growth is 

substantially independent of density ratio because of two 
comj>ensating factors : 

(a) The amplification rate is lower for a higher density 
and given velocity difference. 

(b) The velocity difference is greater for a higher density 
ratio. 

4. The (listurl)ances achieve their terminal amplification 
at a value of A 2 A 1 roughly inversely proi)ortional to the 
frequency of the disturbance. Thus, higher frequencies 
achieve their amplification near the flameholder; lower 
frequencies, remote from the flameholder. 

5. The magnitude of the terminal amplification depends 

upon the inverse of flame slope or roughly The 

amplification that can be a(‘hieved for the maximum of 
2ai^i = () is about (1.1) For example, for 
amplification is eightfold. For Vt/U„=\ /\00, amplification 
is 14,()00-fold. 

6. The disturbance propagates with a phase velocity up 
to GO percent greater than the cold-gas velocity. 

7. The flameholder inhibits the growth of antisymmetric 
disturbances having higher frecpiencies but permits sym- 
metric disturbances to be amplified near the flameholder. 
At 8 diameters downstream of the flameholder, the stability 
picture is substantially that obtained by neglecting the 
flameholder. 


EXPERIMENT: EXAMINATION OF DISTURBANCE 
GROWTH IN A V-FLAME 

INTRODUCTION 

The efl’ect of imposing controlled disturbances on a V-flame 
anchored in a duct was investigated experimentally for two 
reasons : 

(1) To gain quantitative information on the mechanism 
of the breakdown of the flame-generated shear region 
that would explain tlie occurrence of fiame-amplified 
turbulence 

(2) To find a type of excitation that promotes flame 
])ro])agation without causing the sort of shell failure 
that combustor resonance causes 

The method of excitation consists in exciting a column of 
gas in the flameholder support tube. This method offers a 
control on the amplitude when excitation is acconqdished 
witli a speaker, l)ut could be accomplished by means of an 
edge tone excited by the flowing stream for application in 
an engine. Only the speaker excitation is used here. The 
sound issues from a port at the flameholder and disturbs the 
flame at the flameholder. The two types of i)orts used 
operated as follows: 

(1) Sound issued from a port just upstream of one flame- 
holder edge, causing the fluctuating flow about the flame- 
holder to be quasi-antisymmetric. 

(2) Sound issued from the center of the flameholder, caus- 
ing the fluctuating flow at the flameholder to be symmetric. 

The principal disadvantage in this form of excitation is 
that, for some frequencies, modes of burner cavity resonance 
are excited. An objective in the experimental work, there- 
fore, was to determine whether an observed flame-front dis- 
tortion accomi)anying an im])Osed disturbance was due to 
flame-zone instability or to some form of duct resonance. 

No efl’ort was made to determine effects of variations in 
fuel-air ratio or fuel type. The primary variable investi- 
gated was the frequency of excitation. The only attempts 
at correlating the results are the comparison with the results 
of the preceding stability analysis. The chief frustration 
encountered is the impossibility of obtaining a laminar V- 
flanie at velocities where values of VtIUq are sufficiently low 
for the preceding theory to be valid. The reason for this 



Exhaust duct 


Test section 


Calming section 


Figure 12. — Combustor. 
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Air 



is the instability of the boundary layer, which 
separates from the flameholder and undergoes transition to 
turbulent motion quite close to the point of separation. All 
the results are obtained with the flow already somewhat 
turbulent. 

The principal instrument used in measuring flame speeds 


is the photomultiplier probe, first applied in a rigorous way 
to the study of turbulent flame speeds by Clark and Bittker 
(ref. 32). An independent confirmation of the method is 
obtained by comparing combustion efficiency measured by 
the photomultiplier probe and by momentum pressure drop. 

APPARATUS 

BURNER TEST SECTION 

A photograph of the burner (fig. 12) shows the important 
components. Air and propane were filtered, metered, and 
mixed before entering the two-stage calming section (fig. 13). 
In this section the mixture passed through six banks of 200- 
mesh screen and then through a %-inch-radius nozzle plate 
into the ji- by 2-inch (or alternatively 1- by 2-in.) rectan- 
gular glass-walled test section. The top and bottom walls 
of this test section were by J^-inch (alternatively ji- by 
1-in.) brass stock. Tlie side walls were }i-inch Vycor plates 
resting on brackets and held against the }^-inch top and 
bottom pieces by clamps. 

The exhaust system consisted of a O-inch-diameter pipe 
connected to a hood and a i)lower to force the gases out. 
A 4- to 6-inch spacing was maintained l)etween the test 
section and the exhaust duct, and between the exhaust duct 
and the hood; there was no indication of resonance due to 
the exhauster. 



type of disturbance 


Antisymmetric 


Symmetric 


Antisymmetric 


Figure 14. — Flameholders. 
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C- 3833 5 

Figure 15. — Photomultiplier probe. 

Tlie llameholder was a flat liollow paddle 36 to 38 inclies 
long with a 60° gutter on one end and a %-ineh tube leading 
to a piston-type loudspeaker on the other. Tlie three flame- 
holders used in this program are shown in figure 14. The 
flameholders used to supply the disturbances were resonators 
that had resonant frequencies every 190 to 200 cycles when 
coupled to the speaker. The flameholder designed to pro- 
duce antisymmetric disturbances allowed the sound to issue 
from a slot just upstream of one edge of the flameholder; 
the one designed to give symmetric disturbances allowed the 
sound to issue from the center. The flameholdei’s pictured 
were approximately 0.25 to 0.30 inch thick. In one test a 
0.50-inch gutter was tacked onto the 0. 306-inch antisym- 
metric flameholder to examine the effect of flameholder size. 

In no case was the distribution of v' downstream of the 
flameholder measured. It is not known how much the 
experimental distribution diffem from the measured theo- 
retical distribution. The essential features in both experi- 
mental and theoretical disturbances are that antisymmetric 
disturbances cause the flow to flap back and forth at the 
centerline of the duct, while symmetric disturbances cause 
it to dilate about the centerline. 

AUXILIARY EQUIPMENT 

Sound-producing system. — The speaker was supplied a 
sinusoidal signal from a signal generator through a power 
amplifier. The speaker and flameholder could be positioned 
as desired. 

Flow metering. — Air and propane were metered through 
rotameters. 

Hot-wire equipment. — An NACA constant-temperature, 
frequency-compensated, hot-wire anemometer was used 
(ref. 33). Spectra and u' were measured at various locations 
in the duct. Measurements of u' were made primarily with 
a vacuum-tube voltmeter. Two spot checks were made 
against a root-mean-square computer, and the correction 
factor was found to be approximately unity. Spectra were 
measured with a wave analyzer. 


Pressure measurements. — Velocity surveys were made 
using 0.018-inch hypodermic tubing, or 0.030-inch stainless 
tubing, for impact tubes. Measurements were read on a 
water micromanometer. Velocity measurements were made 
at the plane of exhaust with room pressure as static. 

Photomultiplier probe. -The photomultiplier tube was in- 
stalled in a probe as shown in figure 15. Light reached the 
photomultiplier tube after passing through two vertical slots, 
which were and Ke inch wide and 2 feet apart, and through 
filters when desired. Photomultipliers were supplied bv a 
ripple-free regulated d-c source. 

Alternating-current and direct-current voltmetei*s having 
identical input resistances were used where indicated. 

Shadowgraph system. -The optical arrangement used in 
taking shadowgraph pictures is shown in figure 16. The 
light source could be used continuously or as a 4-microsecond 
flash. The camera was focused 1.5 feet in front of the test 
section. The shutter and flash were manually synchronized. 

Stroboscope.— A stroboscope consisting of a slotted alu- 
minum disk and a synchronous motor was used to view the 
disturbed flame in the preliminary phases of the investiga- 
tion. The stroboscope can be seen in the lower right corner 
of figure 12. 
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PROPERTIES OF SYSTEM 

Velocity profiles upstream of fiameholder. — Velocity pro- 
files measured 3 and 12 inches downstream of tlie nozzle are 
shown in figure 17. The fiameholder was positioned so that 
each measurement was taken upstream of the gutter. 
Measurements were also taken at the plane of tlie nozzle, and 
no variation greater than % foot per second from 55 feet per 
second was found. These measurements are for a nominal 


flow setting of 50 feet per second. 

Flame photographs and flame-speed measurements were 
made at nominal velocities of 50 and 100 feet per second in 
K- and 1-inch-wide test sections; the profiles of figure 17 
therefore represent the highest percentage of boundary layer 
encountered. 

Velocity profile in wake; no flame. — Velocity profiles in 
the wake of a 0.306-inch-wide, 60® gutter at a nominal 
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(a) Profiles 3 inches downstream of nozzle. 

(b) Profiles 12 inches downstream of nozzle. 

Figure 17. — Time-averaged velocity profiles with cold flow. Nominal velocity, 50 feet per second. 
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velocity of 46 feet per second are sliowii in figure 18. This 
set of i)rofiIes shovv^s a comparison of the spreading of a wake 
with and without an antisymmetric disturbance at 200 
cycles per second. The only amplitude measurement avail- 
able for this disturbance is a root-mean-square u' of 5 feet 
per second, measured )i inch above tlie lip of the flameholder. 

The effect of other frequencies on wake spreading was 
examined b}^ placing total-pressure tubes 3 inches down- 


stream of the flameholder at the center of the wake and 
varying the frequencies. The results are interpreted in the 
following manner. The velocity is read while disturbances 
are being introduced. The distance from the wake at which 
this velocity is attained without the sound (say x) is found. 
An index of the lengtlis required to accomplish ec[ual degrees 
of mixing without and with the disturbance is then the ratio 
x/3. This ratio is plotted against frequency in figure 19. 





Distance from combustor axis, /, in. 


o Without sound 
□ With sound at 200 cps 
(antisymmetric 
disturbance) 



1 

I 



L 

-J 




O ot 
□ at jT 

'= 81/^' 
= 3" 







Figure 18. — ^ Velocity profiles in wake of 0.306- inch flameholder. 
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Tlie fact that the apparent diffusion rate decreased at 
frequencies greater than 1000 cycles per second indicated a 
decrease in form drag for frequencies greater than 1000 
cycles per second. Static measurements at the base of the 
gutter showed no change in base pressure for the liigher- 
frequenc}^ disturbances. Those lower frequencies that in- 
creased mixing also decreased base pressure (i. e., increased 
form drag). This probably means that the velocity is in- 
creased in the wake for frequencies less than 1000 cycles per 
second. There seems little hope, however, of reducing 
flameholder drag by imposing high-frequency disturbances of 
this kind. 

Turbulence intensity and spectra. -The turbulence com- 
ponent u' in tlie direction of flow was measured with 0.0002- 
by 0.08-inch tungsten wires, using the NACA constant- 
temperature anemometer amplifier. The wire was oriented 
parallel to the flameholder tips. No flame was present for 
these measurements. The results are presented in the 
following table: 


Distance 


^ Distance 

Vertical 


Scale of 

from 

.Nominal 

from 

distance 

'rurl)nlence 

turbu- 

nozzle 

velocity, 

toi> wall. 

from 

component. 

lence. 

Hp, 1 

V, ft/.sec 

//. iti. “ 

station, 

?i', ft /sec 

L, in. 

.r, in. 



in. 



0 

fiO 

1 ~ 

1 .; 

K 

1 0. 300 



100 



. 420 


() 

■/) 

% 1 

H 

1.3S 

0.2 



t ts * 

H 

l.r»s 





K 

.HIS 

. la 




H 

. 07S 




i.i , 

K 

1.021 

. 1.5 



% 

K 

3. 3a 

. 17 


UKl 

H 


2.44 



“ At j=0, 1 / measured from centerline: at. T=0, y measured from flamelioldor lip. 

The spectra measured at several stations were compared 
with spectra found representative for isotropic turbulence 
(refs. 33 and 34). A typical plot is sho^vn in figure 20. The 
scale L is estimated as the one whose spectrum most nearly 
corresponds to the measured spectrum (refs. 33 and 34). 
Tliis scale ^ where determined is given in tlie preceding table. 

The apparent transition from laminar to turbulent flow 
in the boundaiy layer on the flameholder surface occurred if 
the flameholder was about 1 inch or more from the test- 
section nozzle; the disturbance velocit}" in the boundary 
was roughly that in the nozzle itself. With the flameholder 
within 1 inch of the nozzle, eddies were shed at 420 cA^cles 
per second when no flame was present. With the flameholder 
1 inch or more from the nozzle, transition occurred upstream 
of the flameholder lip, and no clearly defined eddy-shedding 
f requeue}" was found in the wake. 

With a flame anchored on the flameholder, no edd}" shed- 
ding was observed without external excitation. The flame 
appeared turbulent for flameholder positions greater and 
less than 1 inch from the nozzle. Therefore, the transition 

8 This is the scale that would be used for calculating a turbulent burning velocity from 
existing theory (ref. 8). It is of interest herein because of the scale of disturbance the flame 
amplifies; 7/«l. 



Figure 20. — Turbulence spectra for flow velocity of 50 feet per second; 

3^(/) = (4h/t/)/{ 1 4-[(27r/f^)h]2/2}, measured Y% inch above flame- 
holder lip. 

to turbulence in the separated boundary evidently affects 
the turbulent flame appearance in about the same manner as 
transition in the attached boundary layer. The lean-blowout 
fuel-air ratio was adversely affected when the transition was 
moved from the attached to the free boundary layer. 

In an effort to obtain a laminar flame zone, some explora- 
tory work was done to move the transition point farther 
downstream of the point of separation from the flameholder. 
Small unheated C 3 dindrical rods (as in ref. 3) placed about 
1 inch from the nozzle entrance were used as flameholders. 
It was possible to move the transition point about 1 inch 
downstream of the flameholder (at 50 ft/sec). Lean blowout,, 
however, was moved to a fuel-air ratio so rich that excessive 
heating occurred in the glass walls, which impaired the use 
of the shadowgraph system. In the transient period im- 
mediately after ignition while the walls were still cool, a 
clear picture of the flame could be observed. Studying the 
flame for a length of time was not possible with this arrange- 
ment. Therefore, turbulent transition at or upstream of tlie 
flameholder lip was accepted. 

CALIBRATION OF DISTURBANCE AMPLITUDE WITH SPEAKER INPUT 

Because of the interaction between the standing-wave 
tube supplying the sound to the flameholder and the stand- 
ing waves in the test section, it was necessary to calibrate 
tlie sound input with the flameholder at the position and 
velocity to be used in the burning tests. For most runs the 
flameholder was placed midway between acoustically open 
ends in the test section in order to suppress velocity disturb- 
ances at the flameholder due to all even modes of vibration 
of the test section. 
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The inputs were set in the following manner: The flame- 
holder position and flow velocity to be used in a burning 
run were set. The hot wire was placed 1/8 inch above and 
1/8 inch upstream of the flameholder lip, to center it over 
the port from which the sound issued. Resonant frequencies 
were excited in the flameholder support tul^e. The input 
to the speaker recpiired to give the desired u' at the wire 
was recorded for each frequency. The wire was removed, 
fuel was introduced and ignited, and the burning run was 
made. The distribution of v' resulting from this sort of 
excitation was not explored. Tlie change in the nature of 
interaction between test sections and the standing-wave 
tube with and without combustion was not explored but 
was presumed to be small, since only about 30 percent of the 
combustible mixture was burned in the test section. 

RATIO OF TURBULENT TO LAMINAR FLAME SPEED DETERMINED BY 
PHOTOMULTIPLIER PROBE 

Photomultilflier tubes have been used to sense transient 
area changes in a laminar flame distorted by sound (refs. 
35 and 36) and to study Bunsen flames that were alternately 
laminar and turbulent by comparing the average intensity 
of light emitted from the flame (ref. 32). Turbulent and 
laminar flames having identical mixture i*atios and feed rates 
produced identical quantities of light (ref. 32). Therefore, 
it appears that the photomultiplier tube, if used with cau- 
tion, can become an average-heat-release-rate meter. 

In the following discussion, the local heat-release rate is 
characterized by the ratio of turbulent to laminar flame 
speed VtIVl. This i*atio is defined as the heat released per 
unit time divided by the lieat that a pair ot laminar flames 
alined with the axis of the duct would release per unit length 
per unit time. The relation is illustrated by considering 
the way the (piantities are measured: 



In the test section shown in sketch (k), flame fronts are 
anchored on the top and bottom of the flameholder and 
appear to span from wall to wall. The photomultiplier 
scans a narrow region of the test section. The voltages read 
from the photomultiplier-tube output are recorded at the 
same fuel-air ratio for a turbulent flame and for a laminar 
flame. The flame-front half-angle is measured, and each 


voltage is multiplied by the cosine of the half-angle. The 
ratio of the resulting product for the turbulent flame to that 
for the laminar flame is defined as TV/Tl- 1^^ the region of 
interest, the cosine of the half-angle was approximately 1 for 
the turbulent flames, so that, in effect, the unit length of the 
flame fronts becomes the unit length of the combustor. 

In obtaining the laminar flame, the same fuel and air 
rotameter settings were used as for the turbulent flame, but 
a large fraction of the mixture was bled off upstream of the 
flameholder. Three banks of 200-mesh screen were installed 
just upstream of flic flameholdei*, and the resulting flame 
appeared laminar, fldie flame did not span from wall to 
wall near the flameholder; therefore, the voltage recorded 
for this region was (piite low. Far downstream of the flame- 
holder the flame front became twisted so that it rejiresented 
a flame depth greater than tlie wall-to-wall deptli. In be- 
tween these two regions was one in wluch the flame appeared 
to touch the walls and to have negligible twist. The voltage 
recorded by the pliotomulliplier probe as it scanned this 
region was nearly constant. The extreme voltages recorded 
in this region envelop a band within which the true voltage 
corresponding to a pair of laminar flames was thought to lie. 

In order to pinpoint an exact value for the laminar voltage 
as well as to verify the method (the zone of uncertainty is 
quite small), an auxiliary determination of VtIVl was made. 
Photomultiplier-probe voltages were read at seven equally 
s])aced stations along the combustor from flameholder to 
exhaust nozzle, and an average voltage was calculated. 
Simultaneously, the momentum pressure drop across the 
flame zone was measured. Witli the lielp of vScurlock’s 
results (ref. 1), the fraction of combustible burned F was 
determined from the measured momentum pressure dro]). 
Then an average Vt was calculated from 

,,_ (Length) (2) (T>) (de])th) 

(l\) (depth) (height) 

where length, depth, and heiglU are the dimensions of the 
combustor. Witli this I>, the measured average voltage, 
and the value for laminar flame speed taken from refer- 
ence 37, voltage corresponding to a laminar flame was 
obtained from 

Vl Vt 

where El ami Et are tlie photomultiplier voltages for laminar 
and turbulent flames. This voltage lies within the zone of 
uncertainty found by direct viewing of the laminar flame. 
These results are reported more fully in the section Exami- 
nation of Method. 

The fact that reference 32 was cited does not mean that 
any photomultiplier calibration was obtained from that work. 
Reference 32 established that (for a given combustor at a 
given composition and flow rate) a photomultiplier voltage 
is proportional to the heat released within sight of the ])hoto- 
multiplier tube and is independent of the existence or non- 
existence of turbulence. The agreement of the two methods 
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for determining VtIVl justifies this use of tlie photomultiplier 
probe. 

It was hoped initially to follow the amplitude of a flame 
wrinlde by measuring the a-c component of the photomulti- 
plier signal. The reasoning was that the flame would be 
wrinkled sufficiently for its local instantaneous slope to serve 
as a detectable index of flame-wrinkle amplitude. This 
slope, of course, would appear as an increased voltage at the 
photomultiplier because the length of the flame cutting across 
the field of view would be increased. The data showed an 
amplitude more dependent on the distance from the duel 
exhaust than on the distance from the fiameholder. The 
sound waves introduced as the soiu-ce of excitation at the 
fiameholder reflected from the exhaust of the combustor. 
Therefore, the a-c component read by the probe was largely 
influenced by the time-varying dis])lacement of local flame 
zones due to plane waves of sound. The filter a-c signals 
for two frecpiencies as a function of distance from the flanu'- 
holder are compared in figure 21. One set was obtained 
from the photomultiplier probe with burning, the other set 
from a hot-wire anemometer without burning. 

The photomulti])lier probe was used only to survey for the 



o 



Distance from fiameholder, in. 


(a) Frequency, 570 cps. 

(b) Frequency, 1500 cps. 

Figure 21. — Comparison of hot-wire and photomultiplier measure- 
ments showing effect of plane waves near exhaust of duct. Flow 
velocity, 50 feet per second. 


effect of distm-bances on local YtIVl values on a time-average 
basis. 

RESULTS 

The analytical part of this paper suggests certain beha ior 
that can be examined experimentally. Accordingly, tests 
are included that furnish the necessary data for a comparison 
with the results of the analysis. One of the primary varia- 
bles studied was the effect of the frequency of excitation on 
disturbance growth. The frequency range of interest hap- 
pened to include some of the resonant frequencies of the 
combustion-gas volume as well as the spontaneous eddy- 
shedding frequency of the fiameholder in cold flow. The 
eflect of combustor j*esonance oj* si)ontan(‘ous eddy-shedding 
was also examined. 

Ortain ])recautionary measures were taken in choosing 
geometry and operating conditions. To guard against 
burner resonance, the combustor was kept short, so that 
only about 30 percent of the combustible was consumed. 
Spontaneous burner resonance sometimes occurred when 
this value was exceeded. The conditions were ke])t within 
the Reynolds number reported in reference 1 to be free of 
spontaneous eddy-shedding when a flame is anchored on the 
fiameholder. Although these measures ])recluded the occur- 
rence of spontaneous oscillations due to these two causes, 
they could not prevent the system from favoring certain 
frequencies within the range of interest. These influences, 
then, received special attention. 

Markstein has shown that instabilities of a laminar flame 
front cause the front to assume a cellular structure under 
certain conditions (ref. 4). The occurrence of such cells 
(‘learly could confuse an interpretation of the results. Con- 
sequently, the data were taken solely at conditions where 
the cells do not spontaneously arise. 

EX.AMIN.4TION OF METHOD 

With the fiameholder 6 inches from the open end of the 
duct and 0 inches from the tunnel entrance, runs were made 
at three fuel-au* ratios at 50 feet ])er second over a range of 
disturbance frequencies. The 0.3()6-inch-wide gutter with 
antisymmetric disturbances was used. The change in im- 
pact pressure with and without combustion was measured 
by a single impact tube upstream of the fiameholder. Since 
the difference is of the order of inches of water and the rate 
of mass flow is constant, it is assumed that the difference is 
entirely due to a change in static pressure. The fraction 
burned was calculated from vScurlock’s results (ref. 1) and 
the measured pressure drop. Concurrently, the light in- 
tensity at seven stations along the flame was measured and 
averaged. 

The fraction burned is compared with the average voltage 
(which is proportional to the average intensity) in figure 22 
by plotting both against the frequency of disturbance for 
the three fuel-air ratios examined. The following table 
indicates the amount of scatter in the ratio of fraction burned 
to mean voltage: 
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Fuel-air 

ratio 

FjEr for frequency, /, cps, of— 

Average 

F/Et 

0 

380 

470 

760 

1500 

0. 0444 

0. 137 

0. 131 

0. 134 

0. 137 

0. 131 

0. 134 

.0515 

. 1017 

.1070 

.1065 

. 1063 

.1035 

.105 

.0580 

.0818 

.0817 

.0863 

.0858 

.0852 

.084 


From the relation 

F 2 Ft-X length X depth 

Et~~ Uo X heigh t X de]) th X Er 

(where length, height, and depth are dimensions of the eom- 
biistor), an average value of VtIEt is calculated. From the 
average values in the preceding table for F/Et at each fuel- 
air ratio, the following values of VtIEt were obtained: 


Fuel-air 

ratio 

VtIEt, 

(ft/sec)/ 

volt 

0. 0444 

1. 115 

.0515 

.880 

.0580 

.709 


For these fuel-air ratios and the inlet temperature and pres- 
sure used, a laminar flame speed was obtained from reference 
37. Then the voltage that would correspond to a laminar 
flame was calculated from 


TT 


VtIEt 


The values of IT from reference 37, the value of El calcu- 
lated from the jireceding relation, and the zone of El meas- 
ured directly by viewing a pair of laminar flames in the test 
section are as follows: 


Fuel-air 

ratio 

Vl 

(ref. 37) 

El 

El by direct 
measurement 

VtIEt 

0. 0444 

0. 74 

0. 663 

0. 658 to 0. 682 

.0515 

1. 10 

1.25 

1.15 to 1.39 

.0580 

1.32 

1. 862 



The El taken from values of Vl! {VtIEt) tliat resulted from 
momentum-pressure-drop measurements fall within the limits 
of El by diiTct measm’ement and are used as the calibration 
of the photomultiplier probe. Local values of \\IVl re- 
ported subsequently are defined by the ratio of voltage read 
at that station to this El for the applicable fuel-air ratio. 

COMPARISON OF FLAME SPEED WITH FLAME APPEARANCE 

Sliadowgraphs were taken at the same conditions used in 
the preceding comparison between average voltage and Ap. 
Some of these shadowgi'aphs, along with the distribution of 





(a) Fuel-air ratio, 0.0580. 

(b) Fuel-air ratio, 0.0515. 

(c) Fuel-air ratio, 0.0444. 

Figure 22. — Comparison of fraction burned obtained by momentum- 
pressure-drop measurement with average voltage obtained by 
photomultiplier survey. Flow velocity, 50 feet per second; length, 
6 inches; disturbance velocity, 2.38 feet per second; antisymmetric 
disturbance. 
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^ L foi* eacli, are presented in figure 23. In most of tlieso 
photographs tlie gross distiirliance eaused by the input sound 
can be identified. If a line drawn along the turbulent flame 
front is faired through all but these gross wrinkles, its length 
divided by the lengtli of a line faired along the undisturbed 
surface should be equal to 

A8rf\ _ Area of flame disturlicd at frequency / 

Area of undisturl)ed flame surface 

A similar laitio is availaide from the average* voltages ob- 
tained by the photomultiplier survey in figure 22, defined by 



The (SfilSu)p^ obtained l)y the former method and the 
(Sd/Su)pm obtained from the photomulti])lier readings are 
compared in the following table: 


Fre- 
quency, 
/, cps 

Fuel-air ratio 

0.0444 

0.0515 


(fO 

\ / pm 


(«).. 


0 

1 


1 

1 

380 

. 18 

1.15 

1.09 

1. 11 

570 

1.304 

1. 185 

1. 19 

1. 12 

770 

1.23 

1.23 

1.10 

1. 18 


1. 10 

1.085 

1.09 

1.07 


In general, tlie corn'lation is good. It may be significant 


that, at 570 cps, (Sd/Su)pm is about 10 percent greater than 
(Sd/Su)pg^ At this frequency, it is relatively easy to measure 
the area; that is, experimental scatter should not be as great 
as 10 percent. The evidence suggests the generation of ad- 
ditional turbulence in the flame, other than the growth of 
the low-frequcncy-imposed disturbance. It is regrettalde 
that a laminar stabilized flame could not be obtained in this 
equipment, because this frequency could conceivably cause 
a transition point from laminar to turbulent flow. 

For each plot in figure 23, u' is a constant, so tliat u' l\\ 
decreases witli increasing fuel-air ratio. Also, VtIVl de- 
creases, for the most jiart, with increasing fuel-air ratio and, 
therefore, witli decreasing u'IVl. 

At a given distance a:, the ratio {SJSu)pm defined by 

\\T^ descends with ascending fuel-air ratio. This is in 

[ 1" 77 t l)u 

keeping with the interpretation in the stability analysis. 
The following table gives tlie values of S^/S^ at 570 cps for 
the three fuel-air ratios: 



Fuel -air ratio 

I, in. 

0.0444 

0.0515 

0.0580 


Sd 

u' 

Nd 

U* 

Sd 

u' 


*S’u 

I V.u 

Su 

Vt.u 

Su 

Vt.u 

1 

1.522 

2. 025 

1. 129 

1.220 

1.070 

1.055 

2 

1. 183 

1.437 

1.205 

1. 107 

1. 122 

. 902 

3 

1. 100 

1. 370 

1. 175 

1.08 

1. 119 

.937 

1 

1. 177 

1.358 

1. 151 

1.058 

1. 147 

.913 

5 

1. 197 

1.352 

1. 129 

1.038 

1. 110 

.873 




r 



i 




ir 



0.0444 0.0515 0.0580 0.0444 0.0515 0.0580 0.0444 0.0515 0.0580 0.0444 0.0515 0.0580 
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(a) No disturbance. 


(c) Antisymmetric disturbance at 570 cps. 


(b) Antisymmetric disturbance at 380 cps. 


(e) Antisymmetric disturbance at 1500 cps. 


(d) Antisymmetric disturbance at 760 cps. 


Figure 23.— Comparison of name appearance and Vt/Vl for three fuel-air ratios for 0.306-inch flameholder. Flow velocity at plane of flame 

holder, 50 feet per second. 
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DISTRIBUTION OF LOCAL Vt/Vl VALUES WITH RESPECT TO VARIOUS 
FREQUENCIES 

The ratio of tuiTulont to laminar flame speed VtIVl is 
plotted against distance from the flameholder in figure 24(a), 
and VtIVl is cross-plotted against frequency at several sta- 
tions in figure 24(b). Great care was taken in keeping 
fuel-air ratio, velocity, and plate voltage to the photocell 
constant throughout the run. The 0.306-inch flameholder, 
modified to prevent leakage about the edges, was used. The 
disturbance velocity measured as described earlier gave u' 
of 2.16 feet per second for all frcciueucies. 

In order to compare these data witli theory, the following 
computations were made: 

(1) From the preceding calibration the fraction burned 
was calculated as a function of combustor length for the 
upper envelope of the VtIVl curves for the disturbed flames. 

(2) An idealized h-ilhi against length was obtained from 
this fraction burned using Scurlock^s results (ref. 1). 

(3) Increments of length over which a given frequency 
has a maximum VtIVl were read from figure 24(a). 

(4) These increments were plotted against holhi and com- 
pared with the theoretical plot taken from figure 10. Recall 



(a) Variation of flame-speed ratio with distance from flameholder. 
Test-section depth, inch. 






Frequency of antisymmetric disturbance, /, cps 


(b) Variation of flame-speed ratio with frequency for various distances 
from flameholder. Test-section depth, inches. 


Figure 24. — Effect of antisymmetric-disturbance frecpiency on flame-speed distribution with 0.30fl-inch flameholder. Flow velocity, 50 feet 

per second; fuel-air ratio, 0.0444; disturbance velocity, 2. Hi feet per second. 
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Figure 25. — Distribution of fraction burned and idcalizc'd /? 2 //b for 
data of figure 21(a). 


that Rusulls of stops (1) and (2) are shown in 

figure 25, and steps (3) and (4) in figure 20. 

Similar data taken at the lower amplitude (^^'=1.44 ft/see) 
gave similar results. 

DISTURB ANCE-DISPLACEMENT GROWTH 

In this section the actual disturbance-displacement growth 
is compared with that of qualitative theory (sketch (f)). 

In an effort to produce low disturbance amplitudes, the 
glass walls of the test section were reduced to 6 inches in 
length and the flameholder was moved to % inch down- 
stream of the nozzle entrance. At this position the speaker 
inputs were calibrated with the hot wire as previously dis- 
cussed, and shadowgraphs were taken for several low am- 
plitudes. As mentioned in the section on turbulence inten- 
sity and spectra, the flameholder boundaiy la\'er underwent 



Figure 20. — Froceiency giving maximum increase in flame speed as 
function of hof/it compared with results of stability analysis. 


transition in the separated layer. The photographs of figure 
27(a) show that, at these low intensities, the apparent widtli 
of the flame was not appreciably altered. The mean flame 
width measured from a direct ])hotograph was used to 
determine 

The displacement h was measured as the distance from 
surrounding crests to a given trough in the flame surface. 
Measurements were made on both sides of the flame for each 
photograph. The data for two amplitudes and two fre- 
(luencies obtained in tliis manner are shown in figures 27(b) 
to (d). The daslied curve is tlie (pialitative curve from 
sketcli (f). 

The knee should occur between {IhlliOmax fiiul (/ioMOcr 
discussed in the section Relation of Theoretical Results to 
Measurable (Quantities. The vertical lines in figure 27 are 
drawn through tliese two values. With the general shape 
of the (lualitative curve known, it can be seen that a line 
faired through the data points would exlhlnt a knee in llie 
region bounded by {h 2 llh)max ifnd (A 2 M 1 ) cr- Comparison of 
figures 27(c) and (d) shows that the knee moves to increased 
values of /lo/Zq as frecpiency is reduced. Both the value of 
A 2 /A 1 where the knee occurs and the variation of this value 
with frequency agree with the results of the analysis. 

DISTURBANCE PHASE VELOCITIES 

From the photographs of figure 27(a), values of CrlUi were 
determined, again using the mean flame width to furnish 
/i 2 //q. The quantity UJlIo was then determined from Tsieids 
results (eq. (3)). Wavelengths were measured from trough 
to trough, and the mean position of a pair of troughs located 
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/ = 955 cps 


/ = 545 cps 








/ = 725 cps 


.8 


.6 


.4 


.2 



0 





Fraction of cross section filled by flame zone, t)2/h^ 


Flame shadowgraplis for disturbance velocity of 0.75 toot per second (b) Disturbance velocity, 0.60 foot per second; frequency, 725 cps. 

at various frecpiencies. (c) Disturbance velocity, 0.75 foot per second; frequency, 725 cps. 

(d) Disturbance velocity, 0.75 foot per second; frequency, 545 cps. 

Figure 27. — Displacement measure on flame front for small-amplitude disturbances. Flow velocity, 50 feet per second; fuel-air ratio, 

0.0484; turbulent disturbance velocity, 0.303 toot per second. 
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Figure 28. — Propagation velocities for sinall-amplitude disturbance. 


the applicable value ot The results are shown for two 

frequencies in figure 28. 

The (laslied lines are the values of CrjUi obtained from 
figure 5(1)), which is based on the results and interpretation 
of the stability analysis. The scatter of the experimental 
data is wide, but tlie magnitude and trends of the data re- 
mote from the flameholder are quite good. Values of CrjUi 
less than unity are obtained for values of ^ 2 /^ 1 ^ ho, wliicli 
implies How reversal within this region. 

EFFECT OF TEST-SECTION DEPTH 

A single run was made to determine the effect of the test- 
section depth on the results. The data obtained witli the 
1 -inch -deep combustor at a velocity of 50 feet per second 
and a fuel-air ratio of 0.0493 are shown in figure 29. The 
turbulence at the center of the duct was half that of the 
/^-iiich duct. Combustor depth seems to have some influ- 
ence on the turbulent flame speed but not on the effect of 
imposed disturbances. 

The undisturbed flame speed in the 1-inch duct reached 
almost twice laminar flame s])eed in a longer time than in 
the /^-inch duct. The effect of imposed disturbance on 
amplification as a function of frequency remains the same. 


Figure 30. —Effect of iiiiiplitude on flame-speed di.stribution. Flow 
velocity, 50 feet j)er second; fuel-air ratio, 0.0444. > 



Figure 29. — Flame-speed distribution in 1-inch-deep test section. 
Flow velocity at plane of flameholder, 50 feet per second; fuel-air 
ratio, 0.0493; J^-inch flameholder. 



Distance from flameholder, jr, in. 
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liowcwer. This can he soon from a comparison of figures 29 
and 24(a). The range over wliich a given frequency domi- 
nates, taken from this run, is also plotted in ligure 26. 

EFFECT OF AMPLITUDE 

The effect of amplitude on single frequencies is shown in 
figure 30, where flame s])eeds for two levels of excitation 
(570 and 960 cps) are com])ared. For the lower frequency, 
the increase A(TV.d/^V, «) varies directly as u' . For the 
higher frequency, the curves cross. Thus, the effect of low 
excitation is felt farther along the flame than that of higher 
excitation. This behavior is not data scatter and was 
checked repeatedly. The same sort of thing occurs at 1500 
cps, although not so markedly. The reason for this anom- 
alous behavior is not immediatel,y ap])arent. 

EFFECT OF FLAMEHOI DER SIZE 

Flame appearance and Vt!Vl compared for two flame- 
holder sizes at a velocity of 50 feet per second and a fuel-air 


ratio of 0.0444 in figure 31, The data for a / 2 -inch flame- 
holder and a variety of excitation frequencies (fig. 31(a)) 
show that the 570-cps excitation is again the dominant one,, 
but not so markedly as before. The effect of flameholder 
size can be summarized from figures 31(b) to (d) as follows: 

(1) Flame s])eed is noticeably decreased with the larger 
flameholder. 

(2) The im])act of the 570-cps frequency is not nearly so 
great in the /^-inch flameholder as in the 0.306-inch flame- 
holder. 

The decrease in flame s])eed with increase in flameholder 
size might be explained by the theory that (1) the wake 
region is longer and (2) the growth of the higher frequencies 
is therefore somewhat inhibited by this long, relatively 
stable wake. 

EFFECT OF VELOCITY 

Momentum-pressure-drop and flame-speed readings were 
taken over a range of frequencies, using a 0.306-inch 
flameholder, a velocity of 100 feet per second, and fuel-air 



/= 570 cps /=l 140 cps /=l360cps /= 1500 cps 



s 


m 

w 


m 



§ 

0.306 


V 

0.500 


Flameholder size 



0.306 0.500 

Flameholder size 



(a) 1^-Inch flameholder with anti- (b) Both flameholders without dis- (c) Both flameholders with disturb- (d) lioth flameholders with dis- 
symmetric disturbances. Disturb- turbance. ance at 570 cps. turbaiice at 1500 cps. 

ance velocities, feet per second: 
u't, 2.22; Utj 1.2; and id, 1.87. 


Figure 31. — Comparison of flame appearance and VtIVl for and 0.306-inch flameholders. Flow velocity, 50 feet per second; fuel-air 

ratio, 0.0444. 
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ratio of 0.0493. Tlio rosults arc shown in fiiuiro 32. The 
metering systems drifted slightly at lliis setting, making it 
diflieult to hold fuel flow steady and causing some scatter 
in the determination of 

The most sensitive frequency at 100 feet per second was 
1700 cps, which is about three times the frequency observed 
at 50 feet per second. Cleai'ly tlie effect is not allied with 
a spontaneous eddy-shedding frequency. Tlie fraction 
burned is not so great as at 50 feet per second. Also, the 
flame s])eed increased sharply u]) to the end of the duct. 
These two factors tend to emphasize the higher frecpiencies 
within the limited duct lengtli, and seem to be directly re- 
lat('d to the instability of the flame zone. The single data 
point at 1700 cps is includecl in figure 20. This run was di- 
rected toward finding a frequency other than 570 cps for maxi- 



mum effect, mainly to prevent duct resonance. Shadow- 
graph ])hotographs of the flames slunv the antisymmetric 
disturbance. 

The possible effect of duct resonance on these results 
must be considered. The modes of resonance likely to occur 
in a combustor are the '^organ-pipe” or longitudinal modes, 
the transverse modes, and a single mode coupling the com- 
bustion chamber with the volume upstream of the nozzle. 
In the configuration studied, the longitudinal modes have 
frecpiencies, I’oughly, of 250 or 500 cps ])lus multii)les of 500. 
The transverse mode's have a minimum frequency of 3000 
cps, so they are not considered here. The last type of oscil- 
lation givc's a frequency of the order of 45 to 90 cps and is 
below the range of interc'st. 

A longitudinal mode could affect the results in two dis- 



Figure 32. — Comparison of fraction burned with average voltage, distribution of IV/ Ul, and shadowgraphs with 0.306-inch flame- 
holder. Plow velocity, 100 feet per second; fuel-air ratio, 0.0403; antisymmetric disturbance. 
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tinct ways: (1) If the llamehokler were placed in a standing 
wave system at a jiressnre nocUy the time-varying velocity 
past the flameholder woidd cause large initial disturbances 
that would be symmetric with respect to the duct axis. 
(2) If t!ie flameholder were placed at a pressure loop, there 
woidd be a cou])ling between the combustor and flame- 
holder resonance that woidd be initially compensated for by 
the in])ut calibration already discussed. If the nature of 
duct resonance shifts the resonant frequencies so that this 
compensation is no longer adequate, the input u' is no 
longer constant for all the frequencies. 

Since the photographs of figure 32 show the disturbances 
to be antisymmetric, the first manner of coupling is rejected. 
To explore the second maimer of coupling, the input u' for 
the frecpiencies on either side of the most sensitive ones 
(i. e., 1500 to 1700 cps; see fig. 32) was increased by a factor 
of 3. In no case was it possible to increase the fraction 
burned as much as the 1500- to 1700-cps disturbances in- 
creased it. Thus, the results shown in figure 32 are due to 
the sensitivity of the flame zone to the imposed disturbance 
and are not due to duct resonance. 

SYMMETRIC DISTURBANCES 

A symmetric disturbance was used for tests with a 0.315- 
inch flameholder, velocity of 50 feet per second, and a fuel- 
air ratio of 0.0444. No separate calibration of input sound 
was attempted, the inputs of tlie antisymmeti-ic-distui-bance 
tests being used. The results are shown in figure 33. Com- 
parisons with flame s])eeds obtained with antisymmetric 
disturbances are given. There was a leak of approximately 
Ke inch past one edge, which provided an additional flame 
surface and increased the apparent mean flame speed near 
the flameholder but did not disturb the results at some dis- 
tance from the flameholder. In general, symmetric and 
antisymmetric disturbances gave roughly comparable shifts 


in flame speed near the flameholder. At a distance from the 
flameholder, antisAunmetric disturbances increased the flame 
speed more than symmetric disturbances. These results are 
in accordance with sketch (f). 

DISCUSSION 

Ten questions raised by the analytical and experimental 
phases of the investigation are discussed in this section. 

1. In what respects do the experimental results and the 
interpretations of the stability analysis agree? — The areas 
of agreement are summai’ized in the following table: 


Interpretation of stability analysis 

Investigation of sensitivity of a flame in a 
duct to transverse velocity disturbance 

(1) The flow is stable to symmetric dis- 
ttu'bances, unstable to antisymmetric 
disturbances. A qualification is made 
near a flameholder where the two types 
of disturbances are equally unstable. 

(1) The flame seems about equally un- 
stable to symmetric and antisymmetric 
disturbances near flameholder, is more 
unstable to antisymmetric disturbances 
far from flameholder. 

(2) In moving along the flame at phase 
velocity, a disturbance of given fre- 
quency grows until the flame width be- 
comes about 0.4 of the wavelength; then 
an apparent decay of the disturbance 
begins. 

(2) A disturbance of given frequency first 
grows, then decays. The critical flame 
width where decay starts is slightly 
smaller than that predicted. 

(3) There is a most-amplified frequency 
having a wave number of 5<2a<Ai<7. 

(3) There is a most-amplified frequency 
which, for the lengths investigated, 
corresponds to 2uohi « 12. 

(4) The phase velocity ranges from a ve- 
locity equal to the cold-gas velocity at 
the flameholder up to 1.5 times the cold- 
gas velocity far from the flameholder. 

(4) The phase velocity behaves much as 
predicted except near the flameholder, 
where it is lower than the cold-gas 
velocity. 

(5) The amplification depends on «7 Ur. u 
for a given d^ 2 /d.r. 

(5) The amplification depends on v'IVt.u 
for a given d/» 2 /d.r. 


The details of the system that would define more exactly 
the shortcomings of the stability analysis as well as the 
adequacies await the develojunent of a transient-velocity 



Dfstance from flameholder, x, in. 


o 


Flameholder 



Disturbonce 

None 


Antisymmetric 


Symmetric 



Antisymmetric 



Symmetric 



None 


/=570 cps 


/= 600 cps 


(a) Syininetric disturbances. Disturbance am- 
plitude, unknown ; input same as for figure 24 (a). 


(b) Comparison of symmetric and antisymmetric flameholders with and without 

disturbance. 


Figure 33. — Comparison of VtIVl obtained with flameholders for producing symmetric and 

Flow velocity, 50 feet per second; fuel-air ratio, 0.0444. 


antisymmetric disturbances. 
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inoasuring techniquo tliat can be used in a flame zone. In 
at least one instance there was evidence that a form of transi- 
tion took place (see Comparison of Flame Speed with Flame 
Ap])earance). The stability analysis cannot indicate wlien 
this takes place. 

2. Would the same distribution of displacement and 

with distance from the flameholder obtain if no growth oc- 
curred in the flame zone? This seems unlikely. For ex- 
am])le, if the in])ut disturbance distribution simulates a source 
at the corner of the llameholder, only near the flameholder 
woidd there be a marked transverse velocity com])onent. 
The highest frequencies would show the greatest increase in 
flame speed. 

3. Do eddies shed? -They certainly seem to. In figure 34 
photographs of a lean flame, both shadowgraphs and direct 



Direct photograph of flame through stroboscope-, frequency, /, 600 cps 



Shadowgraph of excited flame; frequency, /, 760 cps 



C-41538 

Eddy trail 

Fku re 34. — Comparison of photographs of excited Hame ^vith eddy 

trail. 


photographs taken through the stroboscope, are compared 
with a smoke-trace photograph of eddies shedding from a rod 
at a Reynolds number of 73 (ref. 30). The ap])earance of 
flame clearly indicates that eddy-shedding is ])reseiit. 

As the fuel-air ratio is increased, the range over which these 
eddies a])])car near the llameholder narrows. Figure 35 
shows a plot of frequencies at which eddies ap])eared against 
fuel flow from lean to rich blowout. These data were ob- 
tained by exciting the flame at fre(|uencies at which the con- 
stant-s])eed strol)osco])e would permit observation; that is, 
approximately every 100 cps. 

Eddies are characterized by lobes of flame of about equal 
size on cither side of the flameholder. These are contrasted 
to mere sinuosity of the flame. At 700 c])s the lobes ap])eared 
symmetric, indicating duct resonance at this frecjuency. 
At 800 cps the lobes were antisymmetric and were observed 
throughout the fuel-air-ratio range. At the lean fuel-air 
ratios, most of the disturbances produced eddy-shedding. 



Figure 35. — Observations of vortex shedding excited by antisymmetric 

disturltances. 


4. Does the turbulent-flame propagation velocity depend 
upon upstream turbulence? It should. Bolz and Burlage 
(ref. 11) and Mickelsen and Ernstein (ref. 12) have shown 
turbident flame s])eeds to be in fair agreement with Scurlock’s 
theory (ref. 8), omitting flame-generating turliulence effects. 
These data were ol)tained in turbulent “soap-bubl)le^^ flames. 
Thus, in addition to the effect of the instabilities mentioned 
here, the incident turbulence slioidd have an effect. 

The most important point, however, is tliat the instability 
merely acts as an amplifier for a relative!}' narrow range of 
frequencies, the output within limits of the linearized theory 
being proportional to the input. Part of the input is, of 
course, the turbulent motion. The other, and in some cases 
the most important part, is the role played by the sound in 
the duct which, upon hitting the flameholder, produces a 
transverse velocity that may be subsequently amplified. In 
that this sound is related to the incident turlndcnce, the 
turlndence has an effect. 

Some data on an effect of incident turliulence were ob- 
tained wliile exploring the properties of the system. Figure 
36 shows a plot of Et against distance from the flame- 
holder X for two flameholder positions and two fuel-air ratios. 
The table in the section on turbulence intensity and spectra 
shows that the turbulence intensity increases from inlet 
nozzle to exhaust and from the wall to the interior of the 
test section. Taking as representative the value of u' meas- 
ured inch above the flameholder lip, the values of the 
turliulence parameter u' j\\ vary over a factor of 2 for the 
two flameholder positions. It is ap])arent that, everywhere 
except near the flameholder, there is little effect of the initial 
u' IV l values measured. The data are insufficient for pin- 
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(a) Fuel-air ratio, 0.0444. 

(b) Fuel-air ratio, 0.0515. 

Figure 36. — Comparison of relative flame speeds at two flameholder 
positions. Flow velocity, 50 feet per second. 


pointing the eaiiso of this lack of effect. Two surmises 
appear warranted: 

(1) The flame was noisier at the 9-inch position than at 
tlie 6-inch position because of flame-duct interaction; u' was 
measured in the absence of a flame. 

(2) The turl)ulence found upstream of the flameholder 
was due to transition in the houndarv layer. The u' in the 
boundary layer was growing for both positions measured. 
It is thought that the growth of this u' is accelerated when 
the boundary layer separates from the flameholder and ([uite 
rapidly achieves a value that would be relatively insensitive 
to u' prior to separation. 

Ad(litional data on the effect of upstream turlndence were 
obtained b}’- putting three banks of 200-inesh screen just 
upstream of the flameholder. Shadowgraphs of the flame 
in this case with and without a 56()-cps disturbance are shown 
in figure 37. Also plotted in this figure is VtI^l measured 
with and without screens. In the absence of imposed dis- 
turbances, VriyL is lower with the screens than without. 
Although the conditions were set for 50 feet per second, leaks 
were found between the mixing section and the flameholder 
for this run so that the velocity was not known. The reduc- 
tion in VtIVl is thought to be primarily due to the redaction 
in u' by the presence of the screens. 

It is of interest to compare the flame-speed data obtained 
in this investigation with data taken in a similar apparatus. 
Figure 38 shows a comparison of VtIVl obtained by photo- 
multiplier tube with VtIVl obtained by wScurlock (ref. 8) by 
measuring a mean flame width. Because Vl was taken as 
a property of the mixture, the values of VtIVl neai’ the 


flameholder are low. The values of VtIVl obtained with 
the relatively high turl)ulence appear higher than those ob- 
tained by Scurlock. The possible displacement of the abso- 
lute magnitudes of VtIVl due to the different methods used 
would not warrant interpreting these data to give a quanti- 
tative dependence on turbulence intensity. 



With screens 
and no 
disturbance 




Figure 37. — Comparison of flame speeds with and without banks 
of 200-mesh screens upstream of flameholder. 



Figure 38. — Comparison of VtI^l taken in }4-inch- and 1-inch-wide 
ducts with Scurlock’s data. 
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5. How does VtIVl obtained from area measurements 
compare with VjIVl taken from photomultiplier data? -The 
shadowgraphs indicate considerable wriidding of the flame 
surface. In addition, it appears that density gradients exist 
within the outer envelope of flame surface, implying heat 
release there. 

Since no single continuous flame surface is ajiparent, meas- 
uring flame area from, these pictures seems impossible. For 
the sake of comparison, however, some relative surface 
extension should be of inteirst to those who would 
formulate a theory of turbulent flame ])ropagation. Such 
measurements were made in the following manner: The nega- 
tives of a group of photographs were used to project an image 
of the flame on a screen. The image of a 6-inch test section 
was feet long. A map measure was used to find the 
length of the outer envelope of the group of flame images, 
ignoring the obvious density gradients within the outer 
envelope. Lengths were measured on the top and bottom 
surface of each flame, averaged, and divided by ])rojected 
length. The resulting quotient is defined as (tV^ 

In tlie following table, the ratios so obtained are compared 
^vith t he avera ge IV/IT obtained from a photomultiplier 
survey 
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0 
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1.35 
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0. 30f) 
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0 

2. 40 

1.05 

1.49 
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In general, the trends indicated by tlie (Vr/lV)p,« are followed 
by the In view of the obvious internal structure 

and three-dimensional nature of the wiinkles, no argument 
for or against the wriidvled flame as an adequate turbulent 
flame model can be introduced. 

6. How can flame speeds be predicted on the basis of 
these results? — The most iiromising method of ])redicting 
flame speeds would be based on an incident and amplified 
spectra, as recommended by Markstein (ref. 4). For ex- 
am])le: If the incident v' spectrum corresponds to a scale of 
turbulence L of 0.15 inch and each frequency is amplified as 
the theory states, the spectra at /i 2 /Ai = 0.5 and 0.25, for 
1// v 2~20 would correspond to the ones shown in figure 39 . 
Between the two spectra would be a continuousl 3 ^ shifting 
spectrum. What is needed, then, is a model of flame propa- 
gation that would allow for a continuously shifting spectrum 
of turbulence. This model should have the advantage of 
allowing for the appearance of discrete frequencies that 
dominate the rest of the spectrum. 

An interesting point comes up in regard to spectra. It is 
seen experiment all.y as well as theoreticallv that verv little 
effect is derived from exciting the flame at a frequencv less 
than (2/3)(f72Ai) or greater than 4f72/q (see figs. 22, 32, | 



Figure 39. — Spectra of incident and flame-amplified turbulence; 
(/) = (4L/ LO / { 1 + [ (27T/ U) LYP } . 


and 7). Therefore, for the present conditions (i. e., velocit}' 
of 50 ft/sec, 2-in. duct) tlio energy in the spectrum of turiiu- 
lence below 200 cps apparent l.v will not increase flame 
spreading. 

A plot of the percentage of energv in the turbulence at all 
frequencies below / against from reference 38 is shown 
in figure 40. For an assumed scale 7>=(l/6j2/q, the fre- 
quency below the frequency (2/3) (f 72/ii) represents about 
35 percent of the total turbulence energv. This amount of 
energv is contained in disturbances with frequencA’ too low 
to be amplified to any usable extent b^^ the flame zone. 
SimilarW from figure 40, the energy in the spectrum below 
the frequencv /=4 772/0 [or /(Z/f7) ~d.66] is 85 percent of 
the total turbulence energv. Thus, about 15 percent of the 
turbulence energy is at frequencies too high to be appreciablv 
amplified by the flame zone. 



Figure 40. — Turbulence energy distribution for exponential correlation 
coefficient (ref. 38). 
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The scale of isotropic turbulence that would put the largest 
percentage of energy in the freiiuency range (2/3)(672Ai)< 
/< 4 tV 2 /?i can he found from tlie extreme of tlie quantity 
(ref. 36) 

tan-' ( 27 t ( 2 ^ 3 ^ f?) 

or 


In this case the amount of energy within the frequency range 
is 55 percent of the total as compared with 50 percent for 
/.= (l/6)2/ii, 40 percent for L=(l/3)2/q, and 45 percent for 
L=(l/20)2/^i. 

7. Can the results concerning frequency be applied to 
multiple fiameholders'? — If n (lameholders are arranged so 
that the flow Reid can he divided into n identical segments, 
vScurlock’s results hold for each segment. Tlie results in the 
theoretical sections of this report for the stability of the Row 
Reid apply to each segment and rely on the condition that v' 
vanishes at the boundary of each segment. Although one of 
several modes of instability may satisfy this condition, it is 
obvious that the foregoing results do not give any information 
about the other modes. 

To illustrate the occurrence of multiple modes, consider a 
proRle examined by Lin (ref. 22 , pp. 219-220) where 

U=A-\-B sin y 

Jdn shows that the wave numbers for neutral stability of the 
different modes available depend on the number of inflection 
points in the flow and can be given by 



The choice of l=\y yields the following profile, which is 
similar to the one studied earlier in this report for a single 
flameholder: 



Sketch (1) 

For this case neutral stability occurs at a\„=al=Tr^ '6, 0. 

If two llameholdcrs arc present, the following profile may 
be approximated by setting l=2\/. 



Sketch (in) 


For this case the wave numbers for neutral stability of the 
four modes obtained are 





0 


Two new modes appear in addition to the two for tlie single 
proRle. 

Similarly, for three flameholders, 


a\i,=a-=T 


V35 V32 -y yV 


tliat is, there appear four roots in addition to the two for 

/ = X,. 

In short, the single mode investigated in this report is ap- 
jdicable to multiple flameholders, but only as one of several 
modes that combine to give the entire instability. 

8 . Can the excitation at the flameholder permit increased 
combustion rates without destroying the burner walls? — Yes. 
In these tests the sound level at the exhaust ])lane was only 
93 decibels with the flame excited. Mhen the flame is 
excited by the destructive duct resonance, the sound level is 
of the order of 200 decibels (ref. 17). 

In order to apply this type of excitation to an engine, the 
control exerted herein (a speaker system) would not be 
needed. Instead, the disturbance can be introduced by 
suitably placed whistles, mounted on the flameholder and 
driven by the flowing stream. The whistle ideally slioidd 
have a number of strong overtones and a fundamental equal 
to the most sensitive frecpiency f/o/ 2 /i. The overtones 
woidd serve to increase the heat-release rate near the flame- 
holder, while the fundamental would increase the heat- 
release rate far from the flameholder. 

9. What would be the sacrifice in stability and pressure 
drop for such an excitation? — Both these penalties require 
further investigation, l)ut from work that has been done the 
sacrifice does not seem great. For example, when the flame- 
holder was ])laced % inch from the nozzle entrance and a 
6 -inch-long test section was installed, lean blowout at 50 
feet per second occurred at a fuel-air ratio of 0.042. At a 
fuel-ar ratio of 0.0453, the disturbance amplitudes required 
to cause blowout were found for several frequencies. The 
results are shown in figure 41. Minimum disturbance 
velocity was reejuired at 780 cps. At this condition the 
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flainoholdor without flame had an eddy-shedding frequeney 
of 420 cps. The Sti'ouhal number is thus 0.206, agreeing 
well with the results of Roshko (ref. 34). The cold-duct 
resonance fundamental was 1000 cps. The reason for this 
result is not known. The main point here, however, is to 
show the effect of disturbances to be used for promoting 
flame propagation upon flameholder stability. At a fuel-air 
ratio of 0.047, it was impossible to induce blowout with 
speaker inputs as much as 10 times those used in the earlier 
tests. The imposed disturbances will probably narrow the 
stability range slightly but not prohibitively, the amount of 
narrowing depending on the frequencies imposed. 

An interesting observation was made in the study of lean 
blowout. As the lean-blowout fuel-air ratio was approached, 
the flame emitted a low-fre(iuency noise at 87 cps. The 
average heat-release rate dro])])ed abru])tly. A survey with 
the ])hotomultiplier ])robe showed the flame to be inter- 
mittently extinguished for all ])Ositions greater than % inch 
downstream of the flameholder. Typical wave forms ob- 
served at several stations downstream of the flameholder are 
shown in the following sketch: 



Random ])hotogra])hs taken at this condition show intervals 
during which no flame was jiresent in part of the combustor 
(fig. 42). 

It is curious that the flame surface facing upstream at the 
bi-eak in the flame is strong and well-defined both in the 
])hotomulti])lier and the shadowgra])h observations. The 
flame surface facing downstream is conversely ill-defined in 
both methods of observation. Near the ])oint of riqiture the 
photonuilti])lier signal is a])])roximately symmetric. No 
explanation for this is apparent. The most significant 
feature of these data is that this form of resonance accom- 
])anied by intermittent blowoff downstream of the sheltered 
zone was the only form of excitation that reduced the heat- 
release rate in the combustor. On this basis it is tenqiting 
to generalize that all oscillations that do not cause flame 
parting must increase the heat-release rate in the combustor. 

For the frequency of maximum amplification (570 cps), 
an increase of as much as 30 ])ercent in heat release or about 
0.2 percent of the kinetic (Uiergy of the flowing stream was 
obtained in the 6-inch-long combustor for an energy flux 


Sound energy 6 

Kinetic energy of stream ~2500 

In this case the noise level at the exhaust was 93 decibels. 

Another way of increasing the heat-release rate that is 
indirectly suggested is to accelerate the flow locally or alter 
the profile to a more unst able one by use of a convergent - 
divergent combustion chamber. One way these can be 
combined with the introduction of a y' at the desired 
frequency is to construct a combustor with a wavy wall. 
Using this method of excitation would require a new experi- 
mental determination of the desired wavelength and ampli- 
tude of the waves in the wall, but would eliminate any 
effect on llame anchoring. 

10. Can the behavior of axisymmetric flames in cylindrical 
ducts be predicted? —First, the steady-state problem, solved 
l)y Scurlock, of the cylindi-ical flame in a cylindrical duct 
yielded a profile as opposed to the _/\_ profile for plane 
flow. If the arc jiart of the curved profile is assumed 
parabolic, the cylindrically symmetric disturbance equation 
can be solved by elementary functions for symmetric dis- 
turbances. This, however, accomplishes nothing, because 
tlie pi'ofile is neuti-ally stable to symmetric distiu'bances. 
The same result is obtained for the _/\_ profile in the plane 
case. For antisymmetric disturbances, the problem is more 
complicated, and no sim])le solution is apparent. 



Fkjure 41. — ^Distribution with freciiiency of disturbance velocities 
required to produce blowout at fixed velocity and fuel-air ratio with 
0.306-inch flameholder placed % inch from nozzle. Flow velocity 
at plane of flameholder, 50 feet per second; duct length, 5.7 inches; 
antisymmetric disturbances. 
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If the profile is affected by a flameholder in such a way 
that it may more nearly be approximated by a J L. profile 
composed of three straight-line segments and two parabolic 
arcs, a range of instability for synimeti*ic disturbances is 
obtained. This ])i*ofde can be compared with the corre- 
sponding pi'ofile in plane flow _/ \_. Such an analysis is 
given in the appendix. A few calculations are given to permit 
comparison with the data of figure 11, for conditions that hi 
approaches infinity and p is constant throughout. Tlie 
results are compared in the following table: 


Cylindrical flow 



Plane flow 


Sketch (o) 


2b 

2ab 

Plane flow 

Cylindri- 
cal flow 

0. 1 

2. 26 

2. 08 

.2 

2. 55 

2. 42 

1.0 

2. 55 

2.54 


Tlie (*ritical wave numbers agree quite closely. At first 
glance, the results for ])lane flow seem to be a fair approxima- 
tion for cylindrical flow. In order to obtain a rough experi- 
mental verification of this supposition, a 1 )^-inch-inside- 
diameter glass tube and its accommodating nozzle were 
installed in place of the rectangular test section. Tlie 
){- by 0. 306-inch antisymmetric flameholder was retained. 
At 50 feet per second and a fuel-air ratio of 0.047, resonance- 
free operation was obtained at a combustor length of 4)^ 
inches or less. With the flameholder set to give this length, 
the disturbance inputs were calilirated as before to give a 
disturbance velocity of 1.8 feet per second. Tlien the 
burning run was made with the settings indicated by calibra- 
tion. The flame speeds were measured by the photomulti- 
plier probe plus a measurement of an average local flame 
diameter. The results are shown in the following table: 


Fre- 
quency, 
/, cps 

VtIVl 

1" From 
flame- 
holder 

2" From 
flame- 
holder 

4" From 
flame- 
holder 

0 

2.27 

3.31 

3.58 

200 


3. 56 

3. 77 

385 

2.59 

3. 56 

3. 95 

580 

2. 68 

3. 56 

3. 95 

770 

2. 82 

3. 65 

4. 25 

980 

2. 89 

3. 70 

3.80 

1170 

2.89 

3.70 

3. 73 

13()0 

2. 96 

3. 65 

3.73 

1740 

3. 09 

3.65 

3.99 



C-4I54I 


Fioure 42. — Spontaneous oscillations sho^ving intermittent flame 
inteiTU})tion. Flow velocity, 50 feet j)er second; frequency, 85 
cps; fuel -air ratio, 0.0444. 
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The flame speed with this arrangement is almost twice 
as higli as with the plane burner. Part of this increase is 
undoubtedly due to the disturbances introduced by the 
corners of the flameliolder. The fraction l)urned is 20 
percent at the position 4 inches from the flameliolder, 
corresponding to Tluweforo, the frequencies in 

4 inches of this duct may be compared with tliose obtained 
in 6 inches or more of the rectangular duct. The values of 
2aoh\ giving the greatest increase in flame sjieed at the 
termini of the respective ducts are, for Uo of 50 feet per 
second. 



P're- 


2a oft 1 for 

Geometry 

qiiency, 
/, cps 

2ft 1 , in. 

(VT/\'L)miT 

Plane 

570 

2 

11.95 

Cylindrical 

770 

1.5 

12. 10 


Of course, these data were obtained in a single test, 
j)erformed to satisfy cuj’iosity. The agrecunent does make it 
seem likely that a carryover of about 1:1 would be feasible 
in going from ])lane to cylindrically symmetric flow. 

SUMMARY OF RESULTS 

THEORY 

The flow field of a flame in a duct is unstable to transverse 
disturbances having wavelengths 2.5 times the local flame 
width. 

The amount a given frequency disturbance will grow vai*ies 
directly as the fre([uency and inversely as the rate of flame 
spreading dA/d,? and is independent of velocity and density 
ratio across the flame front. The frecpiency / where maximum 
velocity-disturbance growth obtains is /= (2.8/tt) 
where Uo is approach velocity and 2/q is duct width. This 
maximum growth occurs by the time the flame nils G4 per- 
cent of the duct and implies an amplification v' jVo= (1 •108 )i/a' 2 j 
where /\2 = d/c>/da*. Frequencies higher or lower than this are 
amplified less. Higher frequencies achieve their terminal 
amplification at smaller values of lulh\. For example, a fre- 
quency of (l.5/7r)(TV2Ai) has an amplification of v'Ivq— 
(1.07)i/^'^ 2 by the time the flame fills 28 ])ercent of the duct. 

By using a crude model of flame propagation, it is possible 
to predict the ({ualitative effect that amjffified disturbances 
of dis(*rete frequencies will have upon turbulent flame speed. 

EXPERIMENT 

The effects of disturbances imposed at the flameliolder 
are in substantial agreement with those predicted by theory. 


In addition, the results show: 

1. hiddies shed from the flameholder with flame when 
suitably excited. 

2. The ratio of turbulent to laminar flame speed varies 
inversely as flameholder size and directly as the ratio of 
flow velocity to flame s])eed or the ratio of disturbance 
velocity to flame s])eed. 

3. Potentially, the most useful theory of turbulent flame 
propagation would be based on an incident and locally 
shifting spectrum of turbulence. 

4. Gains in local heat-release rates can be obtained with- 
out damage to walls or excessive energy expenditure b}^ 
using acoustic disturbances at a flameholder. 

CONCLUDING REMARKS 

In formulating a theory of turbulent flame propagation, 
the form of i]istal)ility that is discussed in this pa])er must 
be taken into account, as must plane-combustion-wave in- 
stability and flameholder boundary-layer instability. At 
jiresent it is not known whether the latter two forms of 
instability can be encouraged to influence turbulent flame 
pi’opagation more than they encourage it spontaneously; 
from the little evidence that is available, it a])pears that 
they cannot. 

The type of instability discussed herein is of additional 
significance because its effect on the ratio of turbulent to 
laminar flame speed can be increased in some important 
applications. Several methods of increasing this ratio are 
as follows: 

(1) Increase the disturbance velocity at the frequency 
that the flame zone can amplify. 

(2) Increase the ampliffcation rate by locally raising the 
mean velocity in the combustor. 

(3) Increase the amplincation rate by locally varying the 
velocity profile. 

A method of enqffoying all three of these improving factors 
would be to build a combustor having a wavy wall. Al- 
though the second and third methods listed will clearly 
cause an increase in pressure drop, the first method substan- 
tially increases the heat-release rate at the cost of a very 
small fraction of the kinetic energy in the approach-stream 
flow and a slight reduction in the operating range of the 
flameholder. 

Lewis Flight Propulsion Labor.vtory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, October 14, 1957 
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APPENDIX 

AXISYMMETRIC DISTURBANCES IN AXISYMMETRIC FLOW WITH DENSITY JUMPS 


This appeiulLx examines the stability of the flow field aris- 
ing from a flame in a duct for axial symmetry. Tlie results 
are used in the body of the text to compare wave numbers for 
neutral stability of comparable plane and axisymmetric flow 
fields. The symbols are defined where used. 


DISTURBANCE EQUATION 
l^illow (ref. 23) defines the cpiantity \l/ as 

The velocity in the ^-direction is 

u=U+u' 

and in the r-direction, 


(Al) 


v=v 


where 


oz r or 


(A2) 


and the continuity equation is satisfied identically. Pillow 
shows that the disturbance equation then becomes 

[(t/-c) (A3) 


where 




(A4) 


d/*^ rdr r^ 

If the study is restricted to cases where and ? 7 " = 


-U' (i. e., parabolic profiles or flat profiles), then eciuation 


(A3) becomes 


or 


with solutions 


(d/-a2)^=0 


dr-^^rdr 


ip~AxIi{ar) +A2/vi(qj/') 


(A5) 


(A6) 


where 7i and are modified Bessel functions of the fii*st and 
second kind. Therefore, 

= {ar)+A,K,{ar)\ (A7) 


and 


r or 








A 1 j^a / 1 {ar) A~I\ («/') J + 

A2[aA:.(ar)+iA^,(a/)]| (A8) 


where the prime denotes differentiation with respect to the 
argument. 

BOUNDARY CONDITIONS 

Consider the following profile and density distribution: 



Sketch (p) 

Let solutions to equation (A5) be denoted by v?i, (p^y pw for 
regions 1, 2, and 3 with densities pi, po, and ps, respectively. 
As in the plane case, tlie requirements are 

(pi = 0 at hi ^ 

(p\ = (P 2 at ho 

<P2=^z^iK ^ (A9) 

and, because of symmetry, 

<^3 = 0 at r= 0 ^ 


Similarly. 




CM = f M 

\dz/2 \(IzA 


(AlOj 


at ho 
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Since, at r=hi, 

'Pi=-4i/|(a/(|) +xlo/Vi(a/ii) = 0 

tlien 

1 _ 1 

'^7v, («//.) 

At r=ho, (p\ = (P2 <><’. "itli equation (All), 


(All) 


i = ai|^ 


(Pi = ai /,(a/l2)- 


I,(ahi) 


/v,(a/l2) J 


Ki(ahi) 

At r=hs, <P 2 = <P 3 or 

<p2= J^ili(ahs) B2K\{oih:i) 'I 


y 


(A 12) 




The following notation is used: 

a, = h{ah.;)~l].^f;l^ 7l.(a/o) 

Ai (a/ll) 


J 


(A 13) 


Since K^(af)^a> as ae^O, /i(a/’)— ^0 as at — >0, and at /■=0, 

¥>3 = 0, 

0.2=0 (A14) 

Since 

hp (bu' jrbu' , dO\ 

a;=“A d(-+^ a/+" '*') 

tlie pressure eonditions yield at hz (from eqs. (A7) and (A8)) 
— Pi |^(— ac+ ?af i) (pi^J 

= ^P‘ 2 1^(— i) ^a<^2+^ —ia<p 2 J (Alo) 


and at /13 

P3 


^ [— 2 (o:r) + i(af 3)] 

=P 2 ^[—i(oic)+i(aUz)] <p-^—i(od<p 2 (A16) 


BOUNDARY-VALUE PROBLEM 


P^=aK[{ah2) +Y /Vi(a/l2) 




«2 

/3| = /, (a/o) 

0-2 = K\ (ah 2) 

Pi = I {(ahi) 

P4 = Ki{ah3) 

^ 5 =a/I (aW +T" 7 i (ah.,) 
hi 


^6 


=/,(.w{‘!t). 


and 


Ti — (<^^^3) — 

72=a/;(a/i3)+r ^.(a/i3) 
ih 

Setting f "1 = 0 (this merely alters the real part of c by f"i), the 
determinantal equation for .Ij, 7>,, .bb, and (\ is 

«1 —02 

0 0s 0, 

— COt2 C0-^ + 0^ C0-! + 0^ 


0 

-Ti 

0 


d (f 3 (') 0 % ^11 (^ 3 (‘) 0 \Q 0 \\ (f 3 — (‘)^2 


= 0 (A 18) 


One immediate result of equation (A18) is to show tliat, 
if dfVdr vanished at either lu or /13 in region 2, c is real for 
all a's. This can be seen by setting 

06^0S= 




and factoring c from the third row of equation (A18), or 
setting 


0 ,= ^ (aI\)(ahs)+~ lAahd 
Pz P3 


jSio=~ {ocK'd (ahi) +1- Kiiahj) 

P3 h-3 

^n=^/.(«A3)(‘;;„\ 

^'^=P3 A4(a/l3.) 

^n = 0.2=(^)^ = 


(A17) 


0 


and factoi'ing Uz=c from the fourth row of equation (A18). 
h]ither operation leaves a determinantal equation in which 
c is whollv real. 

If 

^(ir 


/(If \ _ 

V dr A2 •' 


(A 19) 


and 


P2 

P3 


if 201 — 06 ^ 
(1202 = 08 

(Jz 0 z= 0 n ^ 


(h04 — 012 
P3 


(A20) 
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then equation (A18) can be reduced to 


COtl — f/2«l 




/3i+c/35 


(C/3 — c)72+— ihyi 


P3 


7i 




C0i2 


«1 

{U3—c)rz+— ihii 

04+ {U3 — C)0l( 


7i 


= 0 (A21) 


Tliis is in the form 


jc+k 771 c+n 
0C+5; sc ^-^e 


= 0 


or 


jsc^— 771 oc^ + c [ ks — 77\ q— OTi) kw — iiq= 0 

and the discriminant of c is 


D + kh"^ + -\ — 2 jwks —2 m qo7i— 

2 jw 771 (j [—2 j70on— 2ksm g—2ks07i-\-A77i okw-\-Ajsnq^ (A22) 

If 1)<CS^, then c is complex and an unstable mode exists. 
The symbols of equation (A22) are defined 


k=^ig-i 

<Xl 

n=0zgz 

o=03--09 

7i 

<Z ~ ^ 3^9 ~ (jA ^ 3 /^3 

P3 7l 

7i 

P3 7l J 


A(23) 


REFERENCES 


1. Sciirlock, A. C.: Flame Stabilization and Propagation in High- 

Velocity Gas Streams. Meteor Rep. 19, Fuels Res. Lab., 
M. I. T., May 1948. ^Contract NOrcl 9661.) 

2. Haddock, Gordon W. : Flame-Blowolf Studies of Cylindrical Flame- 

holders in Channeled Flow. Prog. Rep. 3-24, Jet Prop. Lab., 
C. I. T., May 14, 1951. 

3. Zukoski, Edward E., and Marble, Frank E. : The Role of Wake 

Transition in the Process of Flame Stabilization on Bluff Bodies. 
Combustion Res. and Rev., Butterworths Sci. Pub., 1955, pp. 
167-180. 

4. Karlovitz, B61a: A Turbulent Flame Theory Derived from Experi- 

ments. AGARD Selected Combustion Problems, Butterworths 
Sci. Pub. (London), 1954, pp. 248-262. (Disciussion by K. 
Wohl, pp. 268-274, and by G. H. Markstein, pp. 263-265.) 

5. Berl, W. G., Rice, J. L., and Ro.sen, P.: Flames in Turbulent 

Streams. Jet Prop., vol. 25, no. 7, July 1955, pp. 341-346. 

6. Damkohler, Gerhard: The effect of Turbulence on the Flame 

Velocity in Gas Mixtures. NACA TM 1112, 1947. 

7. Shelkin, K. L: On Combustion in a Turbulent Flow. NACA TM 

1110, 1947. 

8. Scurlock, Arch C., and Grover, John H.: Experimental Studies on 

Turbulent Flames. AGARD Selected Combustion Problems, 
Butterworths Sci. Pub. (London), 1954, pp. 215-247. 

9. Leason, D. B.: Turbulence and Flame Propagation in Premix(‘d 

Gases. Fuel, vol. XXX, no. 10, Oct. 1951, pp. 233-238; discus- 
sion, pp. 238-239. 

10. Mickelsen, William R.: The Propagation of a Free Flame Through 

a Turbulent Gas Stream. M. S. Thesis, Case Inst. Tech., 1953. 

11. Bolz, Ray E., and Burlage, Henry, Jr.: The Influence of Turbu- 

lence on Flame Propagation Rates. Jet Prop., vol. 25, no. 6, 
June 1955, pp. 265-275: con. p. 283. 

12. Mickelsen, William R., and Ernstein, Norman PI: Propagation of 

a Free Flame in a Turbulent Gas Stream. NACA Rep. 1286, 
1956. (Supersedes NACA TN 3456.) 

13. Tucker, Maurice: Interaction of a Free Flame Front with a Tur- 

bulence Field. NACA Rep. 1277, 1956. (Supersedes NACA 
TN 3407.) 

14. Summerfield, Martin, Reiter, Sydney II., Kebely, Victor, and 

Mascolo, Richard W.: The Structure and Propagation Mecha- 
nism of Turbulent Flames in High Speed Flow. Jet Prop., vol. 
25, no. 8, Aug. 1955, pp. 377-384. 


15. Rayleigh: The Theory of Sound. Vol. II. Dover Pub., 1945. 

16. P\‘un, J. B., Forney, II. B., and Garmon, R. C.: Burners for 

Super.sonic Ram-Jets. Ind. and Eng. Chem., vol. 43, no. 7, 
July 1951, pp. 1663-1671. 

17. Blackshear, Perry L., Rayle, Warren D., and Tower, Leonard K.: 

Study of Screeching Combustion in a 6-Inch Simulated After- 
burner. NACA TN 3567, 1955. 

18. Truman, John C., and Newton, Roger T.: Why Do High-Thrust 

Phigines Screech? Aviation Age, vol. 23, no. 5, May 1955, pp. 
136-143. 

19. Kaskan, W. E., and Noreen, A. E.: High-Frequency Oscillations 

of a Flame Held by a Bluff Body. Trans. ASME, vol. 77, no. 
6, Aug. 1955, pp. 885-891; discussion, pp. 891-895. 

20. Tsien, II. S. : Influence of Flame Front on the Flow Fkdd. Jour. 

Appl. Mech., vol. 18, no. 2, June 1951, pp. 188-194. 

21. Fabri, J., Siestrunck, R., and Four6, C.: On the Aerodynamic 

Field of Stabilized Flames. Fourth Symposium (International) 
on Combustion, The Williams <fe Wilkins Co. (Baltimore), 1953, 
pp. 443-450. 

22. Lin, C. C.: On the Stability of Two-Dimensional Parallel Flows, 

pt. I. Quarterly Appl. Math., vol. Ill, no. 2, July 1945, pp. 
117-142; pt. II, vol. Ill, no. 3, Oct. 1945, pp. 218-234; pt. Ill, 
vol. Ill, no. 4, Jan. 1946, pp. 277-301. 

23. Pillow, A. F.: A Review of Hydrodynamic Stability and Its Bear- 

ing on Transition to Turbulent Flow in the Boundary Layer. 
Rep. A. 35, Div. Aero., Council Sci. & Ind. Res. (Australia), 
May 1945. 

24. Schubauer, G. B., and Skramstad, H. K.: Laminar-Boundary- 

Layer Oscillations and Transition on a Flat Plate. NACA Rep. 
909, 1948. (Supersedes NACA WR W-8.) 

25. Pretsch, J.; The Stability of Two-Dimensional Laminar Flow' with 

Pressure Drop and Pressure Rise. Joint Intelligence Objectives 
Agency (Wa.sh., D. C.), July 23, 1946. 

26. Lessen, Martin: On Stability of Free Laminar Boundary Layer 

Betw'een Parallel Streams. NACA Rep. 979, 1950. (Super- 
sedes NACA TN 1929.) 

27. Heisenberg, Werner: On Stability and Turbulence of Fluid Flow's. 

NACA TM 1291, 1951. 

28. Savic, P.: On Acoustically Effective Vortex Motion in Gaseous 

Jets. Phil. Mag., ser. 7, vol. 32, Sept. 1941, pp. 245-252. 


46 


REPORT 1360 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


20. Goldstein, S.: On The Stability of Superposed Streams of Fluids 
of Different Densities. Proc. Roy. Soc. (London), ser. A, vol. 
132, Aug. 1, 1031, pp. 524-548. 

30. Goldstein, Sydney, ed. : Modern Developments in Fluid Dynamics. 

Vols. I and II. Clarendon Press (Oxford), 1038. 

31. Mestre, Andre: Etude des Limites de Stabilite en Relation avec 

la Resistance fles Obstacles a rEcoulement. Combustion Res. 
and Rev., Butterworths Sci. Pub., 1055, pp. 72-86. 

32. Clark, Thomas P., and Bittker, David A.: A study of the Radia- 

tion from Laminar and Turbulent Open Propane-Air Flames as a 
Function of Flame Area, Pkiuivalence Ratio, and Fuel Flow 
Rate. NACA RM E54F20, 1054. 

33. Mickelsen, William R., and Laurence, James C.: Measurement and 

Analysis of Turbulent Flow Containing Periodic Flow Fluctua- 
tions. NACA RM E53F10, 1053. 


34. Roshko, Anatol: On the D(‘velopment of Turbulent Wakes from 

Vortex Streets. NACA Rep. 1101, 1054. (Supersedes NACA 
TN 2013.) 

35. Blackshear, Perry L., Jr.: Driving Standing Waves by Heat Addi- 

tion. NACA TN 2772, 1052. (See also Fourth Symposium 
(International) on Combustion, The Williams & Wilkins Co., 
1053, pp. 553-566.) 

36. Kaskan, W. E.: An Investigation of Vibrating Flames. Fourth 

Symposium (International) on Combustion, The Williams & 
Wilkins Co., 1053, pp, 575-501. 

37. Dugger, Gordon L.: Effect of Initial Mixture Temperature on 

Flame Speed of Methane-Air, Propane-Air, and Ethylene-Air 
Mixtures. NACA Rep. 1061, 1952. (Supersedes NACA TN’s 
2170 and 2374.) 

38. Dryden, Hugh L. : A Review of the Statistical Theory of Turbu- 

lence. Quarterly Appl. Math., vol. I, no. 1, .\pr. 1043, pp. 7-42. 


U. S. GOVERNMENT PRINTING OFFICE: 1958 


